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FOREWORD 
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,y Langlev Research Center, Hampton, Virginia. The project 
lechnlcal Monitor at langley Research Center was Robert E. Shanks, 
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a Morello, Charles E. Knox, and Kathy Sects of Lanploy Resoarc 
Center and Heins Eraberger of S.1SA Ames Research Center arc 
gratefully acknowledged. 
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G0JERATION OF OPTIMUM VERTICAL PROFILES FOR AN 


ADVANCED FLIGHT MANAGEMENT SYSTEM 


John A. Sorensen 
Mark H. Waters 

Analytical Mecnanlcs Associates, Inc. 
Mountain View, California 94043 

SUMMARY 


The objective of this project Is to develop and evaluate one or more 
algor ithn ! and flight management concepts for the on-board minimization 
of fuel or direct operating costs. These concepts are to be used for 
steering a CTOL aircraft in the vertical plane between fixed origin and 
destination airports along a given horizontal path. 

In this report, algorithms for generating minimum fuel or minimum cost 
vertical profiles are derived and examined. The option for fixing the time 
of flight is included in th.e concepts developed. These algorithms form 
the basis for the design of an advanced on-board flight mar agement system. 
The variations in the optimum vertical profiles (resulting from these con- 
cepts) due to variations in wind, takeoff mass, and range-to-destination 
are next presented. Fuel savings due to optimum climb, free cruise alti- 
tude, ard absorbing delays enroute are examined. Finally, the results 
are summarized, and recommendations are made for further work. 

Five appendices are Included which give technical details of optimum 
trajectory design, steering requirements for following these trajectories, 
modeling the aircraft, and off-line computer programs for testing the con- 
cepts. The two computer programs developed are called OPTIM and TR/\GEN, 
and they are available from the Computer Software Management Information 
Center, Barrow Hall, University of Georgia, Athens, Georgia 30601. 
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CoMitrcliU. Airline fllgibr profilAA Art prAtaarly swMtirArAi uainf 
pAAu4o~opeiBlAAtinn pronAdnrAA. Given a reAinnAble tmowlenge o£ the jet 
■ trMia, eoetputer progrem bAced on dyoAtSic progreanlnt teehai<!iuAA teat 
A eeirlee of borlsontAl pethe between the depArture And deetiaAtlon polnte 
to find the one providing nlnimim eoet. The horlAontel pethe teeted 
ere predetermined, end they typlcelly coneiet of a eerlei of Intereon- 
nected eegmente between VOEL/DMB And TACAH etatione. The vertical path 
then coneiete of a fixed Mach number/indicated airspeed profile for cllsib 
and descent combined with a series of iu^tp changes in altitude during 
cruise. The altitude step changes ur>' .'ovemed by predetermJLned fuel 
bumoff rates and air traffic control (*.rc) considerations. This report 
is focused on various aspects of computing Improved vertical profiles. 

This computational procedure should be capable of minimizing fuel 
usage, flight time, or a combination of these variables (direct operating 
cost - DOC). Furthermore, it should be computationally efficient and 
Implemented so that it can be used by the pilot while airborne to change 
the flight profile in case of change in weather conditions, final destina- 
tion, or desired arrival time. Finally, this procedure shoixLd serve as 
a means of automatically driving the aircraft control surfaces and engine 
throttle settings. One element of the work summarized in this report is 
the development of a procedure that can be mechanized on-board to achieve 
the near-optimum path. 

Before continuing with specifics, it is important to mention some 
relevant background concerning the control of flight time. Jet transport 
aircraft and their associated flight schedules have traditionally been 
designed to move passengers and cargo rapidly, and cruise speeds have been 
kept near the maximum design values. Because of air transportation's 
efficiency, this Industry has enjoyed a rapid growth even though airport 
and terminal area congestion is becoming a critical problem. Terminal 
area congestion can produce significant time delays and requirements for 
aircraft holding patterns to absorb these delays. Thus, we have the con- 
dition where aircraft fly at high speeds to terminal areas where they 
then may have to hold for over 30 minutes or even divert to an alternate 
airport because of congestion. In 1975, average arrival delay was from 
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2.9 - 4-6 minuCct at 91 of tha bualaat U.S. airports [5]. It is projactad 
that with tha prasant systam, this dalay will grow to an avaraga of 8.4 
mlnutas by 1990. Tnasa dalays can aasily consuma any fual savad from flying 
a minlmim fuel profila to the terminal area [6]. 

Consider also that to reduce operating costs, the air transport 
industry has been gradually changing flight path characteristics from ones 
that emphasise speed (fast rarvlca; minimum time - dependent costs) to 
ones that tend to minimise a combination of tlme-of-f light and fual. As 
fuel prices continue to climb, it is expected that the nominal flight 
path characteristics will eventually approach those of minimum fuel trajec- 
tories (which take a longer time). A host of flight management avionic 
systems are now appearing which may enable automatic or manual guidance 
along reduced DOC flight paths (e.g., [7]). 

Another element which is affecting the evolution of our alritrans- 
portatlon system is the technological advancement taking place in communica- 
tions capability and digital computer equipment on the ground. The ATC 
management is taking advantage of these advancements to upgrade their system 
through furtlier uses of automation [8], This includes Che eventual develop- 
ment of algorithms and software to command flight profiles which (a) mini- 
mize fuel, (b) provide orderly metered traffic in and out of terminal 
areas, and (c) are conflict free. Thus, the ATC system will be able to 
rmticipate congestion and inform the pilot of what his delay will be be- 
fore ho reaches the terminal area (i.e., the controller will assign the 
pilot a desired t ime-of-arr ival that is conflict free). 

A second important element of this effort has been the extension of 
the optimization procedure to Include the fixed tlme-of-arri’'al constraint, 
fills extension is based on three assumptions; 

a) that aircraft will soon be nominally flown along minimum fuel 
vertical flight oaths, 

h) that delays of variable length will occur at destination ter- 
minal areas because of congestion, and 

c) the AT(.’ system '..ill bt able to predict delay, inform the pilot 
of the magnitude of this delay while he is enroute, and give 
hin: credit for absorbing the delay enroute. 

Hi is studv evaluates the benefit of absorbing delavs enrouti . 
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PrtYloti* tliiiiflcEiit work in ofriwai ¥*rcle«l flight fwtht 

that ninisiiss fual or dlr»:t nitrating eocts ort aiaHMtrlt»l in Infa. 9 «M 10. 
Snhnlti and tngilafcy ivrolopad nathoMtlenl •olutioni of tho m^iaus fotl- 
fixo4 rango-fisHidi tist prohlam (11). Srvin dnrolopMl tho eoneaft of itra- 
togie control [3, 12) in which control of arrival tim Ci.n ha ua^ to radnea 
( )ceaBB ftml that rasulta from dalaya. Also, tha fravioua tima baai^ oMtar- 
ing work of tha MMA Langlay ataff for flight taating |>roflla daBcanti on 
tha Tarminal Ckm figured Vahicla (7CV) aircraft dawmstratad tha practi^lity 
of thaaa aolutiona [13]. Thaaa tasta varifiad pravioua aiamlation studiaa 
conductad by MSA and PAA perBonnal [14). 

Tha NASA TCV t.lrcraft providat a vital link in taating tha faasibility 
and verifying tha axpaetad parforaanca of both tha advanced airboma flight 
management and ATC automation aystama. Tha objective of thla project has 
been to provide both profile optimisation concepts and working models to 
be a part of an advanced flight managamant system design. This design may 
eventually ba flight tested on the TCV aircraft. 

This report is organised as follows: 

a) . In Chapter II, algorithms of methods for generating optimum 

vertical profiles are derived. These mathodB have many options 
including tha ability to fix the time of flight (l.e., the 
tlme-of-arrival) . Tlnese procedures also form the basis for 
specifying on-board mechanization requirements. The algorithms 
have been coded into a computer program called OPTIM. 

b) . In Chapter III, a presentation is made of an evaluation of the 

optimum vertical profiles using OPTIM and a program called 
TRAGEN, TRAGEN simulates an aircraft following a given ver- 
tical path (optimum or otherwise). Sensitivity results are pre- 
sented of the variation in the optimum profile as a function of 
variations in wind, takeoff mass, and range-to-deetindtion. 
Advantages of using free cruise altitude and fixed time-of- 
arrival options are also presented. 

c) . Chapter IV summarizes the study. It also makes recommendations 

regarding further work required to upgrade the computer programs 
and to mechanize the optimization techniques as part of an advanced 
flight management system. 
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App«idic«. A through E ar« pr«.«ntod to glvo tochnicl detail, of trajectory 
optiniaetion, modeling detail, of turbojet engine performance, medium 
range twin>jet tranaport aircraft dynamic, end climb fuel bum e.timatlon, 

and .tearing requirement, for following a reference vertid flight pro- 
file. 
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GQMmATia car thi ommw totigm. ptcariu 

Th« pTlnary vubjMit addr*ts«d k*r« 1« bow to fly tbt siiroraft In 
tb« vnrtloal plana aftar tba borlaontal path haa haan That 

is, tha origin, daatlnatlon, horlaontal routa. atartlng and anding 
altitudas and alrapaada. and takaoff maaa ara flxad, Tha problant than 
la to chooaa tha vartlcal flight prof 11a (altituda. alrapai^, thruat. and 
flight path angla tlnta hlstoriaa) to mlnlmlaa cost. 

In this chaptar. an algorithm la first darlvad which aarvaa aa tha 
basis for ganaratlng an optimum vartlcal proflla. This profUa can althar 
ba ganeratad bafora flight or as a background cc»aputatlon whlla In flight. 
This derivation la basad on trajectory optimisation principles which ara 
outlined In the first part of Appendix A. 

After tha method of generating optimum profiles Is explained, 
example climb and descent profiles are presented for two classes of air- 
craft. Then, modifications to the optimisation algorithms are described 
which allow the user to fix the time length of Che flight. 


Derivation of the Optimization Algorithm 

A set of nonlinear differential equations describing a point mass 
model of the aircraft longitudinal motion Is adequate to describe the 
vertical motion of the aircraft for optimization purposes [16]. This model 
is, for no wind, 

s 

mV •* T cos a - D - mg sin Y, 

mV Y * T sin cx + L - mg cos Y» 

& 

h » V sin Y , (1) 

. a 

X » V cos Y , 

. f 

m “ -w . 


H«r«, the state variables are the true airspeed flight path angle Y> 
altltiKle h, horisontal range x» and suisa a. The control variables are the 
thrust T and angle-of-attack a which are both anplitude linlted. The vari- 
ables lift L and drag D are functions of h, V^, and a. The fuel flow w 
is a function of h, V^, and T. These equations are analogous to Eqs. (A.l) 
in Appendix A. The cost function to be mininixed is of the form, 


J 



( 2 ) 


where and are the unit costs of fuel ($/kg($/lb)) and time ($/hr). At 
this point, the final time t^ is free. The problem is to choose the sequence 
of controls (T and a) that drive Eqs. (1), satisfy the initial and final con- 
straints, and minimise the cost of the flight, as governed by Eq. (2). 


For general trajectory shaping, the flight path angle dynamics can 
be neglected, and it can be assumed that the angle-of-attack a is small. 
Then, the equations of motion of the aircraft in the longitudinal plane 
are 


- (T-D)/m - g sin Y » 

h - V sin Y . (3) 

9 

e 

X “V cos > + V , 
a w 

where the longitudinal component of the wind velocity is included. 

An early attempt at obtaining simplified solutions to the minimum 
t Ime-to-cl imb and minimum fuel-to-climb problems was made by Rutowski [15] 
using graphical methods. He also suggested that the aircraft could adequa- 
tely be represented in such a problem by its specific energy state, 

E - h + V^^/2g . (4) 

9 

The use of this state assumes that kinetic and potential energy can be 
instantaneously interchanged. 
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Bryson, Dssai, and Hoffsian 116) uaad tha anargy stata approxiaaticm 
and dlf farantiatad Eq. (4) to obtain 

E - V^(T - D)/«g (5) 

as the single aquation of sutlon. This assusiad that lift L balanced 
weight tag, and true airspaad Vn was treated as a control variable. 

Now, based on the previous researchers' experience, the specific 
energy and its derivative are defined according to Eqs. (4) and (li). The 
first and second of Bqn. (3) arc replaced with Eq. (5). 

Consider again the cost function of Eq. (2). What is actually de- 
sired is to minimise the cost in flying from origin p^^ to destination P 2 , 
where the time of flight t^ is considered for now to be free. The cost 
function can be changed to have range as the independent variable by 
dividing Eq. (2) by the expression for range rate, or 


5 /•' 

J COB Y + J 

Pi 


2 ^ 
dx 


dx 


( 6 ) 


Pi Hi 

The remaining state equation (Eq. (5)) also is changed as 


dE ^ dE y ^ 
dx dt dt * 

(T-D) / 

mg \V cos Y + V 

a W 


(7) 


This neglects fuel flow as a function of distance. Practical experience 
[17] has shown that as long as the current value of aircraft weight (mg) 
is used in Eq. (7), thi.s simpl It i cat ton is appropriate. 


Equations (h) and (7) are comblnef , using the procedures outlined 
in Appendix A to form the Hamiltonian, 
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H - 


( 8 ) 


dc . , dE 

dJ ^ di 




Cj w + + XV^(T-D)/»f 

V co« Y + V 

a w 


As can ba saan, thla Hamiltonian allows for panallsing of tlma, It 
assumas no particular form for tha fuel flow w, and It takas Into 
account tha horisontal componant of tha wind V^. 

Hote, again howavar, that nona of tha tanas in Eq. (8) ara axpllclt 
functions of tha Indapandant variabla x. Than, according to tha arguments 
assoclatad with Eq. (A. 12) In Appendix A, this Hamiltonian is a constant. 
This greatly facilitates obtaining tha optimtm profile solution. 


Again, assume that the thrust T is a function of the angina pressure 
ratio (EPR) setting 7t. The two control variables can then be considorad 
to be 7t and true airspeed V for this simplified dynamic model. (These 

3 

arguments can be extended so that the conventional controls - throttle and 
elevator deflection - are used). Both ■»! and V have upper and lower 

3 

bounds. By following Eq. (A. 11), the Hamiltonian H ia minimized according 
to 


3H 

9tt 


6n > 0 , 


IM 




0 , 


for all admissible values of tr and V . For an assumed constant value for 

3 

H, Eqs. (8) and (9) constitute three equations with three unknowns - 
TT, V , and X at every point along the optimum trajectory, 

3 


Using the arguments of Refs. J.7-20, assume that the vertical 
flight profile consists of three phases: 

(a) a climb portion where both altitude and energy are 
monctonlcally Increasing with range; 

(b) an essentially equilibrium cruise portion where thrust 
equals drag and lift equals weight; and 

(c) a descent portion where both altitude and energy are 
monotonically decreasing as range-to-go decreases. 


(9a) 

(9b) 
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For th« cruloa porcloo, th« Haaiiltooiaa can than ba writcan aa. 


VLLft 

+ V 

a V 


( 10 ) 


For A glvan cruiaa altltuda and maact Cha EPR aaccing and alrapaad ara 
found by mmerical aaarch tachnlquaa auch that tba valua of thla axpraa- 
alon la ainlnlxad. Thla valua of tha Hanlltonlan la than tha conatant 
valua uaad to clinb to and daacand from thla altltuda. It alao rapraaanta 
tha alnltmm coat par unit dlatanca travalad along tha cruiaa path. 

(Thla Impllas that tha cruiaa altltuda and putaa ara known. Thaaa laauaa 
are alao dlacuaaad later.) 


The conatant valua of can ba aubatltutad Into Eq. (8), and Eq. 
(8 ) can be aolvad for tha negative valua of tha unknown coatate: 


C w + C - H (V coa Y + V„) 
r tea w 

V^(T-D)/mg 


( 11 ) 


Thus, to minimize tha coat function expreaaed by Eq. (8) requires that 
valuea of -n and V be chosen such that Eq. (11) la minimized at every point 
along the climb and descent path. Equation (11) is the key algorithm to 
both offline and on-board computation of the optimum trajectory that mini- 
mlzeit direct operating cost. 


The work of Erzberger et al (17-19, 21) has been based on an alternate 
approach to t' I., problem. He used specific energy E as the independent vari- 
able with the assumption that it always Increases monotonically during climb 
and that it decreases monotonically during descent. Thus, the ranalning state 
variable is the last of Eqs. (3), and Eq. (8) is replaced b> 


H 


dc , , dx 

dE ^ ^ ' 

Cf w + + A (V^ cosY + V^) 

V^(T-D)/mg 


( 12 ) 


11 


Equation (12) cm he co«par«d to Bq. (11) , and It can ba sac a that 
tlMi Haailtonian H and tha costata varlabla X hava ravarsad rolat. Thus, 
thaaa two solutions ara duals of ona anotbar. For anargy aa tha stata 
varlabla, Ersbargar shows, using transvarsallty conditions, that X Is con- 
stant and Is aqual to tha nagatlva valua of tha crulsa cost. In this 
caaa, tha Hanlltonlan varlas in swgnltuda bacausa varlablas within Bq. 

(12) vary with tha suignltuda of tha spaclflc anargy (tha indapandant 
varlabla) . 


Ganaration of Noalnal Optinua Trajactorias 

Bacausa Eqs. (11) and (12) ara assantially Idantlcal In character aiul 
concept. It was convanlant to oMka use of tha coaputar prograa previously 
developed by Ersbargar and Laa (21] to procada with generating optiaua pro- 
files. This prograa is based on tha use of Eq. (12). To avoid confusion, 
tha following description refers to ainlalslng the Haalltonlan of Eq. (12) 
with the costate variable X treated as the constant. The computer program 
uses principles described in Ref. 21, and it Is sunaarised in the second part 
of Appendix A. Tbds program (referred to as OPTIM) , which generates an 
optimum profile between origin and destination points, was extended to hava 
additional desired features, and its specific characteristics (from a pro- 
gramming point of view) have been documented in a companion users' guide 
[22] to this report. It is Instructive here, to present in general terms, 
the elements and steps required to generate an optimum trajectory. 

The following quantities must be input to the program for generating 
points on the optimum profile: 

a) The aircraft initial takeoff altitude, airspeed, and mass, 

b) The final desired aircraft altitude and airspeed, 

c) The range and heading lo be followed by the aircraft from 
the origin to the destination points, 

d) The wind velocity (if non-zero) and temperature profiles (if 
other than standard day) as functions of altitude, and 

e) Tlie values of the constants and in the cost function. 
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Th« program mist have accaaa to nuMrical data (aithar ttblaa or poly** 
itiHilala) which produce valuaa of lift and drag coafficiaota and 
thruat T, and fual flow w aa function* of altitude, true airspead, 
taa^aratura, and EPR aatting. Alao, tor the particular aircraft being 
uaad, the progran miat have the following two aata of eeiplrical equationas 

a) . A Method of gueaalng the approxinate initial cruise Maaa 

baaed on initial takeoff naaa and tha valuaa of fual and 
tine coat; and 

b) . A Method of gueeaing the approxinate landing naaa baaed 

on the final cruiae auiaa 

Finally, the progran nuat have an efficient procedure which can be 
uaed to find the mlniaum value of the Haniltonlan (Eq. (12)) at Vv.ch 
energy level by varying one or nore varlablea (here, and V) over their 
adniaalble reglona. For the OPTIM program, the Fibonacci search technique 
ia uaed. (See Ref. 22) 

It muBt be mentioned that for some caaes, there ia no cruiae section 
of fl'"ht, as is explained in Appendix A. In this case, the constant 
value of the costate variable refers to what would be the cruise coat at 
tne highest point on the trajectory where the profile transitions from 
climb to descent. 


The basic steps tnat the optimum profile generation program follows are 

1) . li.e minimum value of cruise cost (A or II in Lq . ilJ)) is 

i 

evaluated for various cruise masses and with altitude varying 
from sea level to celling height. The results are tabulated 
in a "cruise table". A typical example of cruise cost as a 
function of altitude is shown in Fig. A. 7. 

2 ) . Ba.sed on Initial mass and values of and C^, the amount of 

fuel required tor tin- i 1 imb is approximated (See Kq . (A. 38) and 
Fig. A. 8). This produces cruise mass which is used to deter- 
mine the optimum cruist . ist and altitude irom the cruise t.nble. 
The cruise cost ('* iu . (12); - .(K^,) in Appendix A) is required 

for ninimiza' ;on of Eq. (1.). Also, choosing the cruise cost 
from the cruise table also produces the associated maximum 
values of energy E where the climb portion of flight ends. 
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3). The naxt st«p is to minimis* th* Hmmiltonimn (Eq. (12)) at 
•ach point along tha cliji^ trajectory. Thia begin* at the 
co»r tt*d initial energy and atop* with th* precofsputed final 
climb energy. Th* mmierical procedure procedea as follows: 


a) . A step is made in specific energy - E^_j^ + AE), 

where th* step sixe ii> held constant. 

b) . Th* rang* of acceptable values of and tr are com- 

puted for this energy value . 

c) . Over these ranges, th* Hamiltonian (Eq. (12)) is minimised 

^using the Fibcuaccl technique) with Y assumed to be small 
(cos Y *■ 1)« This, then specifies the values of airspeed, 
lift, drag, thrust, altitude, fuel flow, and energy rate 
(Eq. (5)). 

d) . The time required to make this step is computed as 

At - AE/E . (13) 

e) . The flight path angle Is then computed as 

Y * 8in“^ (Ah/At)/V . (14) 

f) . The value of range ia increased by Ax, where 


Ax ■ cos Y + V^(h))At. 

g) . The aircraft mass is decreased by AW, where 
A'a » w At . 


(15) 


( 16 ) 


It is assumed that final climb energy and Initial cruise energy are 
equal. The above procedure is repeated until the initial cruise 
energy level is reached. This produces an accurate measure of the 
fuel burned during climb and the range traveled during climb. 


4). A final approach (or landing) mass is approximated by 

estimating the cruise range, cruise fuel, and other related 
quantities (See Eqs. (A. 42) - (A. 46)). This procedure also pro- 
duces an estimate of the final cruise altitude, energy, and 
cruise cost. 


5). Ihe procedures for generating the descent profile are now 

followed using the same sequence as for Step (3) but by going 
backwards in time and constraining the energy rate to be 
negative. This is continued until the estimated final cruise 
energy level is reached. This produces an accurate estimate of 
the fuel burned during descent and the range traveled during 
descent. 
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6) . Th« cruise distance is now reconputed by subtracting the climb 

and descent ranges from the total required. This cruise distance 
is used to compute the final cruise mass md the final approach 
(landing) mass, 

7) . With the revised final approach mass, the descent trajectory 

Is reconpjted using Step (5). This also provides refinement of 
the cruise distance, the total fuel burned, and time required 
to travel from origin to destination airports. 

The above description is somewhat simplified in that profiles of 
short range flights (e.g., less than 250 n.mi.) require special computa- 
tlons. Figures A. 4 and A. 5 in Appendix A depict the shapes of the trajec- 
tories for short range for two types of control strategies. (Erzberger 
explains why these two distinct shapes occur in Ref. 17). If both V and 
7T are used as controls, the trajectories have the shapes depicted lu 
Fig. A. 5. For this case with short range, there is no cruise segment, 
and an iterative process is required to obtain the constant value of X 
and the maximum energy state used to compute the climb and descent profiles. 

If only airspeed lb used as the control variable, then thrust is 

set at maximum climb value for climb and idle value for descent. For the 
tri-^jet aircraft model used in this study, this could increase the cost of 
the typical short-range profile by about 1%. However, it simplifies 
mechanization and it reduces the time required to compute the profile. 
Trajectories for this case have the characteristics depicted in Fig. A. 4. 
That is, they have cruise segments, but these segments occur at altitude 
and energy levels lower than the optimal value Again, the program 

must iterate on the correct value of X, and Eq. (A. 48) in Appendix A is 

used to compute the corresponding length of the cruise segment. 

The OPTIM program has been constructed and extended so that various 
options can be exercised: 

1) . No wind or arbitrary input wind, 

2) . Optimize with V and tt or V only, 

3) . Arbitrary cost terms and , 

4) . Fixed or free cruise altitude, 

5) . Twin-jet or tri-jet medium range transport aircraft models, 

6) . Fixed or free time of flight (discussed later), 

7) . With or without a speed constraint of 250 kt below 3048 m 

(10,000 ft.). 
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I). g6«x*t im enili* ov takmoft. for •tunrtiat liii 

ttt« erul»c altitod* c«a b« tct to a ftead ov tlio pvo~ 

gran cam conpmto tha bMt altitude aawiclatad with tha givaa 

etulaa ntaaa. 


Iq addition to ganaratlmg optljnum nominal profUaa» OfPTBl can lia uaad to 
axamlxia tha sonaitivitlaa in trip coat dua to vaxiationa in wind, ranga, 
takaoff mass, imd tha coat tarms and C^. Kicmplas of thaaa waria- 
tions ara givan in tha naxt chaptar. 


Tha computations that taka place in OPTIM for off-line generation of 
optimum vertical profilaa ara also tha same as would be used by an on- 
board flight managamant system. For the on-board mechanisation, the flight 
engineer would enter initial and final altitude and airspeed, initial mass, 
range, and heading. Wind and tmtperature profiles and aircraft mass 
estimates would be updated during flight. The seven basic steps listed 
earlier would be recomputed every few minutes to provide a current reference 
trajectory to be followed. 

The basic OPTIM program was based on using a medium range trl-Jet, 
transport aircraft mode ] 4 It is desirable to extend this program to model 
all major classes of transport aircraft so that the program can serve as 
a flight planning tool for the airlines and others concerned with fuel con- 
servation. A first step In this extension was to add the capability of 
generating optimum profiles for a medium range twin-jet transport aircraft. 
Logic was added throughout the program to account for differences in be- 
havior and characteristics of the two aircraft. 

More details on modeling the aircraft engine are presented in Appendix 
B. The modeling of aerodynamic characteristics of the twin-jet model is 
described in Appendix C. Climb mass estimation procedures are summarized 
in Appendix D. 


Example Optimum Trajectories 

Figures 3 and 4 show plots of various state variables as functions of 
time (or time-to-go) for optimum climb and descent of a tri-jet aircraft 
model traveling 200 n.ml. in range. For this example, the indicated air- 
speed was constrained below 250 kt below 3048 m (10,000 ft). Also, optlmlza- 
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tlon was done by varying V only* with it set to the naximuiii value for climb 
and idle value for descent. The cost constants were $0,265 kg ($0. 12/lb) 
for fu 1 and l?600. 00/hr. As can be seen, when 3048 m (10,000 ft) is reached 
the aircraft levels off and accelerates before continuing to climb. The 
opposite is true for the descent profile. Hote also that on both the climb 
or descent profile, Mach number is generally not constant (contrary to pro- 
files that are specified in aircraft handbooks and normally flown). Snuill 
segments of constant Mach number in Figs. 3 and 4 are due to thrust and drag 
data specified in tabular form at given Mach numbers. This is discussed 
further later. 

Figures 5 and 6 show plots of various state variables as functions of 
time (or timo-to-go) for optimum climb and descent of the twin-jet aircraft 
traveling 750 n.mi. in range. For this example, the Indicated airspeed was 
also constrained to be less than 250 kt below 3048 m (10,000 ft) altitude. 
Optimization was done by varying V and tt, so the thrust could deviate from 

A 

the maximum and idle values for climb and descent, respectively. The cost 
constants were $0.33 kg ($0. 15/lb) for fuel and $600. 00/hr. Again, when 
3048 m (lu.OOO ft) altitude Is reached, the aircraft levels off and acceler- 
ates before continuing to climb. 

By comparing Figs. 5 and 6 with Figs. 3 and 4, it is seen that the 
twin-jet aircraft profiles are similar in character to those of the tri-jet. 
One exception Is that the twin-jet optimum climb profile has a relatively 
large portion where Mach number is held constant. (Again, it is expected 
that this is due to the tabular aerodynamic and engine data.) Also, for 
these paths, the flight path angle is rough because of the inherent assump- 
tion in the optimization theory that either altitude, airspeed, or both can 
be abruptly changed to obtain a change in energy. The flight path angle 
can either (a) be smoothed after it is generated to provide a more flyable 
path, or (b) the optimization process can be constrained by allowing only 
small changes in flight path angle between energy steps. Converting the 
aerodynamic and engine input data to polynominal form may also remove this 
problem. 




Figure 3 . 


State Variables for an Optimum Climb Profile of a 
Medium-Range Tri-jet Aircraft on a 200 n.mi. Range Flight. 
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Figure 6. Concluded 
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Figura 7 couparas tha altituda va. raiifa profllaa for thraa caaaa of 
tha trl-Jat aircraft travallng 200 n.ai. For thaaa. both and ir wara 
variad. Howavar, a haod wind and tail wind wara addad in two of tha caaaa 
to produce tha diffarancaa. Tha wind aagnituda profile ia alao ahown in 
Fig. 7, and it ia rapraaantativa of tha Denver wind naaaurad during Auguat, 
1977. (Thia wind ia generally froa the waat.) 

Table 1 nakaa varioua range and coat compariaona of tha thraa pro- 
fllaa ahown in Fig. 7, baaad on valuaa of of $0.138/kg ($0. 0628/lb) and 
of $500. 00/hr. It la aean that tha head wind cauaad a total coat in- 
craaaa of $24.95 (or 4.3X) while the tall wind dacreaaad tha coat $23.59 
(or 4. IX) compared to the no wind optimum caaa. Data, auch aa praaantad 
in Table 1, are printed for each run of OPTIM ao that tha trajectory coata 
and other characteriatica can ba readily evaluated (again, aaa Ref. 22). 



Figure. 7. Comparison of Minimum Direct-Operating-Cost Flight Profiles for 

tlie Tri-jet Aircraft Model. 
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Tabl* 1. CoH|Nirlwm of Flight Costa for 200 q.bI. Esngs with 
Varying Ulnd Cooditlona 




Case 



)(o Wind 

Bead Wind 

Tall Wind 

Pual Usad (kg (lb)) 




CllBb 

1945.8 

1 2013.0 

1805.2 

(4289.60) 

' (4437.82) 

(3979.66) 

Daacant 

272.2 

295.1 

284 . J 

( 600.19) 

( 650.57) 

( 628.13) 

Total 

2218.0 

2308.1 

2090.1 

(4889.79) 

(5088.39) 

(4607.80) 

Range (n.oil.) 




Climb 

115.51 

115.90 

107.36 

Descent 

85.24 

85.78 

94.60 

Total 

200.75 

201.68 

201.96 

Time (min: sec) 




Climb 

17:31 

18:33 

15:41 

Descent 

15:05 

15:33 

16:11 

Total 

32:36 

34:06 

31:53 

Cost ($) 




Climb 

414.51 

432.40 

379.91 

Descent 

163.29 ' 

170.35 

174. 3C 

Total 

577.80 

602.75 

554.21 

Cost/Distance 

($/n.mi.) 


! 


Climb 

3.59 

3.73 

3,54 

Descent 

1.92 

1.99 

1.84 

Total 

2.88 

2.99 

2.74 


Fixed Tlme-of-Arrival 

One of the problems associated with increased use of air travel and 
air transportation of car^o Is the increased congestion that is occuring 
at the major hub airport terminal areas. Aircraft arrive In the terminal 
area on a somewhat random basis. If the aircraft arrive in too great a 
number, in a short period of time, they are currently placed in holding 
patterns by ATC. Then, they are instructed to land so that airport arrivals 
are at an acceptable rate. The holding pattern approach to regulating the 
landing rate has several shortcomlngb : 
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1. Th« holding pattnrn is *xp«nBiv« to the airline* in terns of 
increased fuel requirssients , excess wear on the aircraft, and 
increased crew flying tine. 

2. The holding pattern delays the passenger, possibly causing 
missed connections. 

3. The risk of collision is increased by holding aircraft in the 
already crowded terminal areas. 

4. Maintaining several aircraft in safe holding patterns increases 
the air traffic controller's work load substantially. 

One solution to thlc problem is to regulate the time of arrival of 
each aircraft such that they do not arrive in random groups but rather 
at a systematic scceptbhl.'ie rate. This is the ultimate goal of the FAA's 
metering, spacing, and flow control projects. This requires that the air- 
craft have the means to compute and regulate the flight profile so that 
the desired time-of-arrival is met. 

The airlines also like to fly their aircraft with a scheduled, pre- 
determined trip time. This time is selected to produce acceptable operating 
costs and desired service to the airline customers. Thus, there is a need 
to control the time length of a particular flight. With these motivations, 
OPTIM was modified to include this capability. 

If tlme-of-a lival t^ is fixed, then the objective is to minimize the 
total fuel burned .n urc.vellng from to in time t^. This Xj equivalent 
to minimizing die function 

J = / w dt (17) 

o 

with the additional constraint that the range traveled is fixed. 

The technique used in OPTIM to compute the profile that produces fixed 

range and fixed tlme-of-arrival is to fix ^ and to iterate on the value 

r 

of C^. For each value of chosen, a given tlme-of-arrival t^ results. 
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As is Incrsassd* Chs valus of dacraasas, and vica varsa. Thus, an 
outar loop has baan addad to tha OPTIM logic to adjust to produca tha 
dasirad t^. 

Ona constraint that had to ba built into tha progran is tha ranga 
that tha variabla could span. As incraasas, tha anphasis is placad 
on achiaving minlntsn tlnia. As fual prlcas incraasa howavar, this is not 
practical. On tha othar hand, as dacraases, tha tima-of-arrlval in> 
craasas. This is accooipaniad by a slowing down in cruise speed. A 
practical lower bound on cruise spaed is that value where fuel 

rate w is minimised. If more time is required in cruise to reach the 
destination at soise desired time t^^, then the aircraft should make up 
the difference by path stretching at the cruise speed 

This posed the question of what should the lowest value of be 

such that the resultant optimum cruise speed is This value was 

determined as follows. 


The cost of flight during cruise is expressed by setting the Hamil- 
talnion value in Eq. (10) to the coefficient X, or 


A 


w + C 

■ V +’ V 
a w 


t 


(18) 


Mathematically, the question is what value of will provide a trimmed 
cruise condition such that the trimmed cruise speed will equal 

At optimum cruise speed, the cruise cost A is minimum, or 



0 . 


If this is also the speed where fuel rate is minimum, then 



0 


(19) 


( 20 ) 
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at this spaad. Applying Eq. (19) to Eq. (18) ylalds 


^ V + V Jv ... ^ „ v2 

a a w a (V^ + V^) 


By substituting Eq. (20) into Eq. (21), we gat 

Cf w + Ct - 0, 


( 21 ) 


or 

C. . “ C, w . 

tmin f min 

Note, that at this value, from Eq. (18), 


( 22 ) 


A . « 0. 

rain 


(23) 


Mathematically, this means that there is no cost to cruise at this condition. 


a 

Figure 8 shows the minimum fuel rate w for the tri-jet as a funct ’.on 
of cruise altitude and cruise mass. As can be seen, flying at minimum 
fuel rate causes a large variation in optimum cruise altitude as cruise 
mass is varied. Figure 9 illustrates how optimum altitude and cruise 
airspeed vary as functions of cruise mass for the minimum fuel rate con- 
dition. 


The results of Figs. 8 and 9 are used to compute the value of 
and tlie associated cruise altitude in OPTIM. For example, for the tri- 
jet with a cruise mass of 54.43 tonne (120,000 lb) and a of $. 33/kg 

($. 15/lb), w . is 2.585 tonne/hr (5686 Ib/hr) . This produces a , of 
min tmin 

-$852. 90/hr at an altitude of 9.25 km (30300 ft). The range on for 

this aircraft is -$777.60 to -$1042.50 for cruise mass varying from 50 - 
66 tonne (110,000 lb to 145,000 lb.) 

Figures 10 and 11 illustrate how tlme-of-arrival and fuel used vary as 
varies between +$300/hr and -$900/hr. , for various combinations of initial 
mass and range-to-go for the tri-jet. In these example plots, cruise altitude 
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Figure 8. Minimum Fuel Rate for the Tri-Jet Model as a Function of 

Cruise Altitude and Mass. 


Is fixed at 10 km (33000 ft). The cost of fuel is set at $. 33/kg ($. 15/lb). 
Note that in Fig. 10, the time-of-arrival has a non-parabolic shape; this is 
different than would be expected. This indicates the presence of a conver- 
gence problem in OPTIM which is discussed shortly. 


The fixed i ime-of-arrival option requires repeated generation of new 
cruise conditions associated with each successive value of time cost 
tried. Thus, the OPTIM logic was set up so that this process was as 
efficient as possible. No short range flights were assumed so that 
cruise altitude was always reached. The table representing cruise condi- 
tions (the cruise table), spanned altitudes from 6 km - 12 km (20,000 ft - 
40,000 ft). If a fixed cruise altitude was chosen, only that altitude was 
used in generating the cruise table. 


■jT 
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Figure 9. Optimum Cruise Altitude and Airspeed for a Tri-Jet Aircraft 
Flying at Minimum Fuel Rate as Functions of Cruise Mass. 


The cruise table consists of two parts; 

1. The cruise speed setting which minimizes A, as defined by Eq. (18), 
and 

2. The cruise speed setting which minimizes w. 

The second part of this table is generated only once at the beginning of 
the run because the conditions do not depend on a particular C^ setting. 

The program first sets C^. to zero. (This corresponds to the minimum 
fuel profile.) The resulting tirae-of-arrival T^ is then compared to the 
desired arrival time If T^ is larger than then flight is too 

slow. In this case, C^ must be positive. Next, trial values of C^ of 
$300/hr and $600/hr are used to obtain arrival times of T^^ and T 2 . A 
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Figure 11. Fuel Burned Variations as Functions of Takeoff Weight 
and Cost-of-Time C^. Cruise Altitude Fixed at 10 km 

(33000 ft). 500 n.ml. Range. 


quadratic equation is derived to pass through the points (T^, , T 2 ) 

as a function of C^.. This quadratic equation is then solved for the 
value of that will produce t^^. This process is repeated using the 
last three values of tlme-of-arrival until convergence is reached. 


If the t ime-of-arrival T is less than t^,, then the minimum fuel 

o fd 

profile is too fast. In this case, must be negative. The next trial 
value the program uses is from Eq. (22). This produces the maximum 

time of If t^j is greater than a holding pattern (or path 

stretching) is required to make up the additional time. If t^.^ is less 
than ( * Tj^), then linear interpolation is used to obtain the next 

value of Cj.. This produces t ime-of-arrival T 2 . After that point, three 
values of t ime-of-arrival are available to compute a quadratic equation 
and solve the C^.. This procedure is again repeated until convergence is 
reached . 


As predicted from Fig. 10, a convergence problem did arise in testing 
the fixed t ime-of-arrival option. This is illustrated in Fig. 12 which 
shows the variation in time-of-arrival for a 61.4 tonne (135000 lb) air- 
craft with varying from S300/hr to $600/hr. As can be seen, there is 
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Figure* 12. 


Illustr.ition of Time*-of-Arrival Convorgonco Problem 




about a 100 sac gap in arrival time around of $400/hr. A quadratic 
optimisation algorithm was triad in place of tha original Fibonacci search 
technique to determine if tha optimisation loutine was the source. This 
produced some Improvement but did not remove the problem as is shown in 
Figure 12. Tightening the tolerances on the cruise trim conditions also 
did not remove the problem. 

The problem was traced to the way tha coefficient of drag was 
determined by OPTIM In computing the cruise table. Currently, is 
computed as a polynominal function of lift coefficient C^, with Mach 
number as a parameter, as Illustrated in Fig. 13. Linear interpolation 
is then used to obtain for the given cruise Mach number. 

Figures 14 and 15 show detailed plots of cruise cost A of the tri- 
jet as a function of cruise speed, with cost of time set at $300/hr and 
$600/hr. Cost of fuel la $. 33/kg ($. 15/lb). Cruise altitude is a para- 
meter. Note from these figures that: 

1. There Is a double minimum at the upper altitudes. The OPTIM 
program is based on there being only one minimum A at each 
cruise altitude. The optimization routines used also are based 
on the assumption of a single minimum. 

2. Tlie minimums are along the Mach 0.76 and 0.82 lines. Tlie Mach 
0.80 lino tends to pull the cost curves up creating the double 
minimums. 

The drag coefficient should be a continuous function of Cj and Mach. 

Tlie minimum cruise cost should not always be at either M of 0.76 or 0.82 
as Indicated in the cost curves. The cost curves should be smooth with a 
single minimum at each altitude. Thus, the aerodynamic data, as currently 
stored in OPTIM, produces a convergence problem in attempting to obtain 
fixed t ime-of -arrival with tolerance on the order of 2 sec. 

An effort is now underway to correct the convergence problem. Steps 
being taken are: 
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1. 


Th« Cp v» Cj^ and M data ara baing checkad carefully to anaura 
that thara ara no klnka or benda in tha curvea. 

2, Linear interpolation ia being eliainated. Polynoainal func- 
tlona ara being developed to anaura that all data vary con- 
tinuoualy and anoothly. 

The aama conatrainta need to be applied to modeling tha propulaion data. 
The leaaon learned la that numerical optimization in baaed on the natural 
exlatance of continuoua, amooth aurfacast and any departure from thia type 
of model will produce a faulty reault. 

Deapite the current axiatence of the convergence problem, there ia a 
logical aolution to its removal. Alao, aa aeen in the next chapter, 
having fixed time-of-arrlval capability can produce a significant flight 
fuel savings. 
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PROGRAM VERIFICATIOH AHD SEH8ITIVITY RESULTS 

The previous ch«pe*r and Appendix A describe an efficient -way in 
vhlch near-optlnum flight profiles can be generated wltlu>ut using t late- 
consuming numerical techniques. Instead of iteratively solving the two- 
point boundary value problem, assumptions are made so that the dynamics 
are simplified to two state variables. One state variable (either range 
or energy) becomes the independent variable, and the other is included as 
the single state in the Hamiltonian. Thus, the problem of solving for 
minimum cost profiles reduces to algebraic minimization of the Hamiltonian 
at each point along the profile. 

Now, although this method is a convenient way of generating a trajec- 
tory, the resulting trajectory must be verified by using a more accurate 
model of the aircraft dynamics. Verification implies that; 

1) The reference trajectory that is generated must be flyable 
when the full aircraft equations of motion and constraints 
are taken into account. 

2) The trajectory cost (Eq. (2)) as predicted by OPTIM must be 
essentially identical to that experlen. cd by simulating more 
complete aircraft <.quations of motion. 

Thus, the first objective of this chapter is to describe a companion 
program to OPTIM which was developed for verification of the optimization 
results. This program is referred to as TRv\GKN (for trajectory g ; oration). 

With OPTIM and TRi\GEN as computer tools, the user has the capability 
to study the ,;haracter ist ics of optimum profiles in great detail and to 
examine alternate wavs these profiles can be implemented on-hoard. Irajec- 
tory characteristics are obtained bv exercising OPTIM 's options and by 
making sensitivity studies with OPTIM and TRAGIIN which are the subjects of 
the second part of this chaf ‘r. 
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Optimization Varification 


To achieve optimization verification required the development of a 
companion computer program to C^-TIM which can be uaed to simulate the 
longitudinal trajectory of an aircraft commanded to follow the reference 
path output from OPTIM. Some details of this program (TBAGEN) are pre- 
sented in Appendix E, and a separate user's guide for this program has 
also been written [23] • 


In addition to verification of the optimization program's results, 
the TRAGEN .'rogram has the following utility: 


1) . It provides a means for testing guidance laws for steering 

the aircraft to follow the input reference trajectory. 

2) . It enables study of the effect of following an Incorrect 

reference trajectory. For example, if the OPTIM results were 
based on one particular wind profile and initial aircraft weight, 
and a different weight and wini profile actually existed, the 
TRAGEN simulation would allow assessment of the effect of these 
errors on trajectory cost. 

3) . It can be used to determine the flight cost that would result 

from the aircraft being commanded to follow a reference trajec- 
tory suggested in the manufacturer's aircraft handbook. For 
example, for climb, handbook reference trajectories usually 
consist of following a constant indicated airspeed until a 
given Mach number is reached. Then, the reference trajectory 
follows this fixed Mach number until the reference cruise 
altitude is reached. 

4) . It can be used to test perturbation control schemes for removing 

the effect of wind gusts and other non-nomlnal performance 
sources (navigation errors, transient temperature profiles, 
non-standard engine performance). 

5) . It is expandable to test candidate on-board mechanizations of 

a system cf equations for generating the near-optimum vertical 
profile. 


A five state-variable model of the aircraft is currently used in 
TRAGEN to simulate longitudinal motion. State variables are altitude, 
altitude rate, longitudinal range, airspeed, and aircraft mass. Currently 
neglecued are the rapid transient dynamics of throttle response, angle-of- 
attack, and pitch rate (6^, a, q). The throttle is assumed to be set so 
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that laaxloiuffl thrust is achieved during climb and idle thrust is used 
during descent. The cruise phase is not 'Imulated. The altitude control 
variable is taken to be the angle-of>attack which has maximum and minimum 
limits. This degree of sophistication is adequate for testing the OFTIM 
results. 

The TRAGEN program can readily be expanded to include throttle 
dynamics and short period dynamics. This would be required for further 
study of autothrottle and autopilot design to steer the aircraft to follow 
input reference trajectories. The control variables would be throttle 
position and elevator deflection for this expanded capability. A require- 
ment for implementing this expanded simulation would be to obtain the 
necessary stability and control derivatives to complete the dynamic model. 


The specification of reference profiles used in TRAGEN is based on 
using altitude as the independent variable for climb and range-to-go to the 
destination as the independent variable for descent. For climb, the 
reference trajectory consists of specifying airspeed and flight path angle 
(with respect to the air mass) as piecewise linear functions of altitude. 

At the 3048 m (10000 ft) point, tl»e aircraft is commanded to level off and ac- 
celerate until the airspeed is reached where the climb should again continue. 


To generate the control law to follow the commanded climb profile, a 
linear perturbation model was made of the dynamic equations. 


T cos a - D(a,h,V ) - W sin Y * mV , 
a a 

T sin a + Lfa,h,V ) - W cos y “ n>V Y • 
a a 


(24) 


The perturbation equations and transfer functions from Eqs. (24) are 
given in Eqs. (E.IO) and (E.ll) in App'^ndix E. Here, it is assumed that 
a perturbation .Su to the nominal angle-of-attack can be used to obtain 
the desired perturbations in flight path angle (6 y) and airspeed (>^>Vg) 
maintain the aircraft on the desired reference climb profile. 
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Because both couraanded values of airspeed and flight path angle are 
assutned to vary linearly with altitude, they are ramp functions, and a 
Type 1 control system is suggested. The original climb control law was 
thus set to be 

K K 

- (1=1 + J?) + (Kj + -}) (Y^ - Y) (25) 

where V and are the commanded reference values of and y . Kj^, K 2 * K^, 
and are control gains. Hie commanded values were derived according to 
Eqs. (E.12) and (E.13). 

Experience with this control law revealed two points; 

1) The Integral control gain K 2 on airspeed error was noi needed 

and did not particularly improve performance. Thus, it was 
nominally set to zero. 

2) Constant values of the gains Kj^, K^, and could be selected 

to provide good, stable performance throughout the climb phase. 

Thus, there was no need to have altitude dependent gains pro- 
grammed . 


No attempt was made to select the control gains so that the perturba- 
tion response was optimized. The main objective was to obtain a set of 
gains which caused the aircraft simulation to track the input reference 
trajectory with only a small amount of error (wliieh was accomplished). 

(Jain selection for control response optimization should remain as a task 
to be conducted when an actual autopllot/autothrottle is being implemented 
cind the complete aircraft dynamics are being considered. 

Evalua^on of Optimum Profi les Several optimum trajectories were 
computed by using OPTIM and then subsequently used as inputs to drive the 
TRA(;EN simulation. I'or example, Table 2 presents a comparison of OPTIM 
and TRAGF.N results at the end of the climb portion of a tri-jet aircraft 
model liaving an initial mass of bl23b kg (135000 lb), and traveling a total 
range of 150, 225, and 275 n.mi. with no wind. As can be seen, the match is 
exceptional. The same degree of comparison was found using the two pro- 
grams for different range and different initial mass flights. 


A4 


Table 2. Comparison of Optimisation and Trajectory Generation 
Program Climb Results for Different Range Plights. No Wind. 
61236 kg (135000 lb) Takeoff Mass. 


Rongtt 

Fu«l Burn - 

Tlwt 

- 8«c 

AlttCuda 

■ (f"t) 

Alrspaad 

min 

■(ft/.) 

RMIg* 

m 

" (ft) 

n. Hi. 

WTIH 

TRAGEN 

[ OPTin TRACEir 1 

TJFnS 

TuACkh 

Opt IM 

TRAGEfT" 

OPTIH 

“TRXtrra 

150 

1405 

1405 

682 

680 

7597 

7597 

230 

229 

125127 

125054 

(3097) 

(3097) 

(24924) 

(24926) 

(756) 

(752) 

(410522) 

(410281) 

22S 

1931 

1931 

1060 

1060 

9760 

9760 

244 

244 

211441 

212936 

(425B) 

(4258) 1 

(32021) 

(32021) 

(801) 

(800) 

(693704) 

(698609) 

275 

2075 

2073 

1185 

1178 

10319 

10321 

244 

244 

242116 

241033 

(4574) 

(4570) 
1 

(33855) 

(33861) 

(802) 

(800) j 

(794345) 

(790790) 


The results with head and tail winds were not as close. Table 3 pre- 
sents a comparison of OPTIM and TRAGEN results at the end of climb when the 
wind profile of Fig. 7 was used as both a head and tail wind for the 225 n.ml. 
range flight. In both cases, the more detailed simulation from TRAGEN shows 
that it takes a longer time period (14-22 sec) and greater range (4000 - 5000 m 
(14000 - 15000 ft)) than predicted by OPTIM (However, range traveled is not 
significant for climb.) The biggest discrepancy is the l.bX extra fuel re- 
quired for the tail wind case. In the future, the modeling simplifications 
of OPTIM and the steering accuracy of the TRAGF.N control law should be in- 
vestigated to resolve this point. 


Table 3. Comparison of Optimization and Trajectory Generation 
Program Climb Results in the Presence of Winds. 225 n.mi. 
Range. 61236 kg (135000 lb) Takeoff Mass. 


Type 

Fuel Burn - 

Time 

■■ Sec 

Altitude “(f*) 

Airspeed 

m/a 

■(ft/s) 

Range 

n 

■ (ft) 

Wind 

OPTIM 

TRAGEN 

OPTIM 

IKAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

Head 

2084 

2082 

1180 

1202 

10294 

10362 

244 

244 

225177 

229818 

(^594) 

(4591) 

(33772) 

(3399b) 

(802) 

(802) 

(738770) 

(753997) 

Tall 

1701 

174.7 

897 

911 

9145 

9074 

233 

235 

183238 

187433 

(3751) 

(3851) 1 

(30003) 

(29770) 

(766) 

(772) 

(601174) 

(614936) 
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For the descent portion of the flight, the more likely variable - range- 
to-go to the destination point - was used as the Independent variable for 
specifying the reference path. In this case, the descent altitude was 
specified as a function of range. This defined the inertial flight path 
angle that the aircraft should be on which Is defined by Eq. (E.16). In 
the current version of TRAGEN, no attempt was made to control airspeed 
during the descent. The control law to compute the perturbation to the 
angle-of-attack was of the form 

6a - (K3 + (Yj^ - Yj) , (26) 

where Yj^ Is the commanded inertial flight path angle. 

Despite the lack of control of true airspeed during descent, the descent 
steering provides the same degree of matching performance between OPTIM and 
TRAGEN results as was experienced during the climb simalation (where true air- 
speed was controlled) above 3048 m (10,000 ft). Table 4 compares OFTIM 
and TRAGEN results at about 3048 m (10,000 ft) for the tri-jet descending with 
idle thrust from 10668 m (35,000 ft). The three cases are with no wind and 
the Denver head and tail winds shown in Fig. 7. The biggest error is in 
airspeed where variations of +5.8 to -1.5 m/s (+19 to -5 ft/sec) are seen. 


I'able 4. Comparison of Optimization and Trajectory Generation 
Program Results Descending from 10668 m (35000 ft) to 3048 m 
(10000 ft). 


Type 

Fuel 

Burn - 

Tl*e-tO' 

-go- • 

Altitude 

IB 

■ (ft) 

Alrepeed 

m/s 

"ift/s) 

Range-to 

m 

■ (ft) 

Wind 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

OPTIM 

TRAGEN 

No Wind 

117 

117 

444 

447 

3154 

3148 

149 

151 

56704 

56598 

(259) 

(258) 

(10348) 

(10328) 

(490) 

(497) 

(186036) 

(185688) 

Head 

110 

110 

444 

449 

3154 

3148 

149 

148 

55465 

55494 

(242) 

(242) 

(10348) 

(10328) 

(490) 

(485) 

(181973) 

(182066) 

Tall 

123 

124 

440 

441 

3154 

314B 

149 

155 

57911 

57856 

(271) 

(273) 

(10348) 

(10329) 

(490) 

(509) 

(189998) 

(189816) 
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These seme cases arc conpared again at the altitude of 1.8 n (6 ft) in 
Table 5. The head wind case has a -3.4 si/s (-11 ft/sec) difference in true 
airspeed Which causes 7 sec tlsie and 4.5 (10 lb) fuel burn differences at 

this point. The tall wind case has a 4-5.8 m/s (4-19 ft/sec) difference in 
airspeed which causes - 14 sec time and -1.8 kg (-4 lb) fuel bum differences. 
These differences could possibly be lessened with true airspeed control added. 


Table 5. Comparison of Optimisation and Trajectory Generation 
Program Results. Descending from 10668 m (35000 ft) to 1.8 ro 
(6 ft). 


typ* 

li"" 

Pu«l Bum - 

1 

O 

BBI 

tSSBBSa 



IS38S&B!SH| 

wind 

IbmM 

1 


EiZHOi 

KiniaB 




OPTIM 



■ni 




l.S 

-3.1 

124 

126 


-49 



msSM 



(*) 

(-10) 

(407) 

(413) 

u 

(-157) 


230 

234 



l.S 

-5.2 

124 

121 


-96 


(506) 

(516) 



(6) 

(-17) 

(407) 

(396) 

u 

(-314) 


263 

241 



1.8 

-4.3 

124 

130 


-85 


(535) 

(531) 



(8) 

(-14) 

(407) 

(426) 


(-280) 


Despite the differences between OPTIM and TRAGEN results, It Is believed 
that the match between them Is very good. That Is, the results provided by 
OPTIM (fuel burn, time expired, trajectory followed) can be concluded 
to bt accurate, and sensitivity runs based on using both programs will also 
produce accurate conclusions. A better match between the two programs can 
be obtained by minor adjustments to both programs. 

E valuation of Handboo k P rofiles It is also desired to use TRAGEN 
to evaluate and compare the cost of climb and descent profiles as specified 
In the pilot handbook with those generated by OPTIM. Thus, it was required 
to simulate flight along profiles specified in the pilot handbook. This 
required adding the capability of computing a typical handbook reference 
trajectory as part of the TRAGEN code. 

Typical reference climb trajectories consist of the following sequence: 
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1. Climbing at 2r>0 kt indicated airapeed (VIASl) to 3048 m 
(10,000 ft.) (250 kt below 10,000 ft la an ATC maximum air- 
apeed conatraint). 

2. Leveling off at lOUOO ft and accelerating to a new climb 
indicated airapeed (V1AS2) (e.g., 320 KIAS). 

3. Climbing at VIAS2 until a Mach number (M3) (e.g., M ” 0.78) 
is reached. 

4. Continue climbing at Mach number M3 until cruise altitude is 
reached . 


During this climb, full throttle is typically used. 

The handbook descent profile is the reverse of the above sequence with 
throttle set to idle. Thus, during descent M3, VIAS2, VIASl are used with 
the same constraint below 3048 m (10000 ft). 

To compute the reference profile which follows that specified by a 
handbook, the following sequence of computations are used In TRAGEN. First, 
the altitude is incremented as 

h - h, + Ah , (27) 

b 

Here, the subscript b refers to the value computed during the previous 
cycle. In TRAGEN, Ah Is currently set to 152.4 m (500 ft). Next, the 
desired value of true airspeed is computed from the specified indicated 


airspeed or Mach number M by 

p =“ fj(h) (28) 

■r - f,(h) (29) 

e i 

y - p/3092. 4Tg (30) 

a = 65.76 y/T (31) 

M « (5.((((V^^g/29.)/p + 1.)^^^ - 1.))^^^ (32) 

Vj, » M a (33) 
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H«r«, p and ara prasaura and tamparatura which are conputad in a aub- 
routina aa functlona of altitude. Than, air danaity p and apaad-of-aotmd 
a are computed aa functiona of p and T^. If Mach nunbar ia given, Bq. 

(32) la akipped, and Inaiaad 

VjAs " 29 .(p((1. + 0.2M^))^-5 - ^ 

ia used to compute the indicated alrapead. 

Next, the apeciflc energy at h la 

E - h + V^2/2g . (35) 

Then, for climb, the EPR is set at the maximum value, and thrust Th and 
fuel flow w are again obtained from subroutines as 

Th - f3(M,h) , (35) 

w - f^(M,h,Th) . ( 37 ) 

The approximation is made that lift L balances the aircraft weight w, or 

L - W . (38) 

Then, the coefficient of lift c; is 

ij 

s L/(l/2 p S) . (39) 

where S is the aircraft reference area. Then, another subroutlnt' is used 
to compute the drag coefficient as 

* ^5 * ( 40 ) 

From this, drag is 

D i 1/2 P S . (41) 

Then, energy rate is computed as 

E i (T - D) V^/W . (42) 
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The Increaent of tine to reach thlK new aitate le 


( 43 ) 


(44) 


AR - (V-, + V ) COH V At/2 . 

The weight is Incremented as 
W “ •” w At . 

The altitude is held constant to accelerate from the Initial speed V 

‘ o 

to VlASl and to accelerate from VlASl to VIAS2. 

For convenience, the descent reference profile is computed backwards 
In time from the final altitude to the cruise altitude. Then, the data 
points are reordered so that the descent reference profile la followed In 
real time. 

Comparison of Optimum and Handbook Reference Profiles 

One value of the TRACKN program Is that it can be used to compare the 
fuel and direct operating costa of various profiles followed In climbing 
to and descending from cruise. This capability gives the user a valid way 
of establishing tlie potential uorth of anv given fllglit management system. 

lahle 0 , taken from Kef. 24, indicates the handbook climb and descent 
schedules for several types of commercial aircraft, lor example, a 727-100 
aircraft climbs (above lOOOO ftl at 140 kt (IAS) until obtaining 0.7H Mach. 
The climb is continued at 0.78 Maeli until cruise altitude is reached, 
t.'ruise is at t'.80 Mach, descent Is also at 0.80 Mach until obtaining 340 kt 
(IAS). Then, descent is continued at 340 kt until 3048 m (10000 ft) is 
reaehed. At that time the aircraft levels off and deeelerates to 250 kt 
before continuing to descend. 


(45) 


(46) 


The flight path angle to reach this state la 

Y - Bln"^ (2(Ah/At)/(V + V ) . 

b 

The change in range between the two points is about 


At i (E - E^^)/E . 
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Table 6. Typical Aircraft Characteristics 
as Specified In Manufacturer's Handbooks. [24] 
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TRAGEN w«« u«ed to co«p«r« fuel end tine cohl < of the handbook type 
profile! with optlauii prof lice generated by OPTIM. Figure 16 Illustrates 
the speed/altitude profiles followed In a comparative example using the 
tri-jet aircraft. Shown are a handbook reference profile, an intermediate 
reference profile, and a miniousi fuel profile generated by OPTIM. The 
intermediate reference profile was selected to have the characteristics 
of the handbook profile but to have speeds more nearly equal to the mini- 
mum fuel profile. 

Table 7 presents the fuel used, time, and range covered in flying the 
three climb and descent segments illustrated in Fig. 16. An incremental 
cruise distance was added to the Intermediate and minimum fuel climb pro- 
files so that the same range would be covered for the climb comparison. 
Similarly, a cruise increment was added to the handbook descent profile so 
the same range would be covered for the descent comparison. (In Table 7, 
the second number after each + sign is the added cruise portion.) 

As can be seen, the minimum fuel climb profile uses 169 kg (371 lb) 
less fuel (7.8%) than the handbook reference climb profile. It also takes 
36 sec less time and covers 15.338 n.mi. less range to reach the cruise 
conditions. Note however, that the minimum fuel profile begins cruise at 
Mach 0.77 rather than 0,78. Also note that if the cruise segment is added 
to the minimum fuel climb so equal range is covered, the improvement in 
fuel is cut to 52 kg (115 lb), or 2.3X. The minimum fuel climb profile is 
essentially equivalent to the intermediate profile in that it only uses 6 kg 
(13 Ibl less fuel. 

The descent profile shows a different result. Jlere, the handbook descent 
profile taken alone uses 27 kg (60 lb) less fuel and 267 sec leas time to 
descend than the "minimum fuel" profile. If the range traveled adjustment 
is made, however, the faster handbook descent requires 57 kg (125 lb) more 
overall fuel. This includes the estimated 84 kg (184 lb) required to main- 
tain cruise at 0.8 Mach for 10 n.mi. fhus, this faster descent consumes about 
20% more fuel, overall, than the rainlroura fuel profile (although about the 
same mass of fuel is gained in using the minimum fuel descent as in using 
the minimum fuel climb.) 
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Table 7. Comparison of Steering Techniques on Cost Performance for 

Tri-Jet Aircraft. 500 n.mi. Range 10 km (33000 ft) Crulae 
Altitude. 



Altitude 


_ Handbook Reference Profile 

Fixed KIAS/Mach Profile 

MininuM Fuel Profile 



Airspeed - kt 


Figure lb. Comparison 


of Optimum and Handboor Reference Type Profiles 


other Profile Coaperleone end Sencitivlty iitudlee 


Opttaua Profile Cherecterletlce In conducting the eeneitivity 
enelyeie of the optlimin profilec, there ere tvra iteae that are of par~ 
ticular liitereat: 

1) . The eatabliehment of how the characteriatice of optinum profiles 

change with variatlona in initial aaas, range-to-go, wind pro- 
files, and other variables affecting the aircraft performance. 
Knowing these variations affects how an on-board algorithm should 
be constructed to account for measurrd changes in the flight 
conditions. 

2) . The analysis of the effect that errors in the assumed flt^tt t con- 

ditions (l.e., different initial mass, different wind profile) 
have on the cost of flying a particular profile. For example, 
if assuming the wrong initial mass of the aircraft has little 
effect on the overall cost, then the algorithm can be constructed 
so that initial mass Is assumed to be a nominal value. 

The purpose of this section is to begin to address these Items for wind, 
initial mass, and range variations. Both OPTIM and TRAGEN are required 
for this purpose. Other variables that may have an effect are temperature, 
lift coefficient, drag coefficient, and thrust variations as well as 
measurement errors of altitude, flight path angle, airspeed (or Mach 
number), and range-to-go. These variations should be addressed at a future 
r Ime. 


it is assumed for climb that Mach number and flight path angle given 
as functions of altitude specify the optimum profile to a given cruise 
altitude. Thus, bv presenting plots of these variables as functions of 
altitude, the effect of the variations to the nominal conditions is directly 
seen . 


Figure 17 shows the variations in the climb profile as the total range 
is varied, in steps of 25 n.mi. from 75 n.mi. to 275 n.ml. for the tri-jet 
model. As can be seen, there is no variation below 3048 m (lOOOO ft). 

Above 3048 m (10000 ft), the flight path angle is initially higher and the 
climb Mach number is initially lower as the total range is increased. But 
it Is seen that these curves are essentially parallel, and they merge into 
a common curve for range exceeding 175 n.mi. Thus, it would be relatively 
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Mach Number Flight Path Angle - deg 



Altitude - km 
a) Flight Path Angle 


—4 

40 

"T2.5 



Altitude - km 

b) Mach Number 


Figure 17, Effect of Range Variation on Optimum Climb Flight Path 
Angle and Mach Nuiiiber 


«£.sy to medal thaaa curves as a polynoraloal function of altitude with para- 
meter variation due to change in range. 

Figure 18 shows the variation In the profile as initial mass is 
varied. Again, there appears to be a parallel offset with a greater effect 
on flight path angle. Tnere is no effect on Mach number below 3048 m 
(10000 ft).) Also, again it can be seen that this variation would be easy 
to account for in specifying an on-board reference trajectory. 

The biggest uncertaitity that would be experienced during the climb 
and descent would be the variations in the longitudinal wind. To study 
this effect, the two wind profiles of Fig. 19 were used as program Inputs. 
The Denver wind is the same as is shown in Fig. 7, but Fig. 19 also has 
the wind headings indicated at discrete points. The "triangle" wind is 
a constant 270° heading wind varying linearly in magnitude with altitude 
at a rate of 1 knot per 305 m (1000 ft). These wind profiles were assumed 
to act as both ta*' and head winds by setting the aircraft heading to be 
90° or 270°. 

Figure 20 shows the effect of each of these wind profiles on the 
optimum climb trajectory compared to the nominal profile having no wind. 

As can be seen, the Mach number variation with altitude is essentially 
a parallel offset. The effect on flight path angle is mainly seen going 
from zero to 3048 m (10000 ft) where a fan out of +0.3° is seen. Above 
3048 m (10000 ft), the variation appears more as a steady offset. Thus, 
again it appears that accounting for various wind profiles can be done 
in the on-board system in a relatively simple way. 

Figure 21 shows the effect of each of the wind profiles on the optimum 
descent from 10 km (33000 ft). There, the plots show altitude and Mach 
number as functions of range-to-go to the landing point. From the first 
plot, it is seen that wind causes a + 10 n.ml. range variation in where 
the optimum descent should begin. This is predictable based er> knowing 
how far the air mast moves due to the wind during the descerx ; riod. Also, 
there is up to a + 0.05 variation in the optimum Mach number at a 'unction 
of range-to-go. From the earlier results presented, it was shown that if 
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Figure 18. Effect of Initial Mass 
Flight Path Angle and Mach Number 


Variation on Optimum Climb 


30 

Wind Speed - kt 

Figure 19. Wind Profilc-e I'sed to Study Sensitivity Effects on Optimum 

Profiles. 
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Figure 21. Effect of Wind Variations on Optimum Descent Altitude and 
Mach Number as Functions of Range-to-go. 


the altitude-range profile was followed closely down to 3 k» (10000 ft), 
the Mach number would probably follow Che speed profile associated with 
Che given altitude range profile. Thus, from an implementation point- 
of-view, the Mach variations seen above 3 km (10000 ft) in Fig. 21 are 
not of concern. Below 3 km (10000 ft), only one profile needs to be stored. 

The results seen in Figs, l?'*^! are limited in scope, in that all possible 
variations are not addressed. For example, the effects of variations in 
cruise mass and cruise altitude for the descent trajectories should be 
explored. However, the one thing that can be concluded is that variations 
of total range, initial mass, and wind profiles have predictable effects 
that can be easily used to modify the characterization of the nominal 
optimum trajectory. Thus, the process of computing the optimum climb-des- 
cent profiles does not have to take place on-board the aircraft. Instead, 
a nominal profile plus modifications to account for off-nominal parameter 
changes can be pre-computed and stored in the aircraft flight management 
system. This provides a very simple on-board method of computing the 
reference profile. 

The second aspect of the sensitivity analysis was to study the effect 
of flying a non-optimal profile. This effort involved using TRAGEN to 
simulate following a given incorrect optimum profile when a different 
profile should have been used. Two error conditions were investigated: 

a) Following an optimum profile specified for the incorrect 
ir.itial mass 

b) Following an optimum profile specified for the incorrect wind 
prof lie. 

Only the climb phase was investigated. 

For initial mass errors, four runs were made as shown in table 8. As 
can be expected, the increase in assumed Initial mass causes a decrease 
in the optimum cruise altitude. Thus, if the aircraft is lighter than the 
assumed amount, it will climb faster to the lower altitude. This results 
in both time and fuel reduction from what was predicted to reach cruise. 

Cases la and Id in Table 8 were of this nature, and they show a 13-14% re- 
duction in fuel and a 14-16% reduction in climb time for an initial error 
of -6804 kg (-15000 lb) in initial mass. 
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Tabls 6. Initial Mass Error Cases Considered ’•Climb Profile 


Run 

Initial Maas 
- kg (lb) 

Assunad Optimum Profile 
Mwa - kg (lb) 

Change in 
Climb Fuel 
kE (lb)-(%) 

Change in 
Climb Time 
Sec - (%) 

la 

61236 

68040 

-274 



(135000) 

(150000) 

(-604) (-13.1%) 

-159 (-14.4%) 

b 

61236 

54432 

+329 



(135000) 

(120000) 

(+726) (+18.7%) 

223 (+22.9%) 


68040 

61236 

371 



(150000) 

(135000) 

(+817) (+19.2%) 

247 (+23.3%) 

d 

54432 

61236 

269 



(120000) 

(135000) 

(-592) (-13.9%) 

-167 (-15.8%) 


On the other hand, if the initial mass is assumed too small, the 
heavier aircraft will attempt to climb to a higher than nominal altitude. 
Cases lb and Ic of Table 8 show this condition for an initial mass error 
of +6804 kg (+15000 lb). The result was a 19% increase in required climb 
fuel and a 23% Increase in required climb time. It was also seen that 
the aircraft did not achieve the required cruise airspeeds when the in- 
correct optimum input altitude was reached. Climb fuel estimation accuracy 
is important for choosing the optimum cruise conditions. 

Based on the results of Table 8, it is seen that the Initial mass is 
an important parameter to be entered into the computations. It is better 
to assume too large an initial mass than vice versa. 

For the climb wind profile errors, four more cases were studied where 
optimum climb profiles based on the triangle wind profile of Fig. 19 were 
used. The results are shown in Table 9. As can be seen, if the tail wind 
is greater than assumed (Cases 2a and 2d), more time and possibly more fuel 
are required to achieve cruise conditions. Likewise, if the head wind is 
greater than assumed (Cases 2b and 2c), less time and possibly less fuel 
are required, although t'le results are mixed. It is seen that neglecting the 
wind profile can vary fuel cost up to + 1.5%. 
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Tabl« 9. Wind Profila Error Casas Consldarad - Cllisb Profile 


Run 

Actual Wind 
Profile 

Aaaumed Profile 
for Computing 
Optimum 

Change In 
Climb Fuel 
. (lb) - (X) _ 

Change in 
Climb Time 
Sec - (%) 

2a 

No Wind 

Head 

-1 (-2) (-.05%) 

+7 (+0.7%) 

b 

No Wind 

Tall 

+11 (+24) (+0.6%) 

-3 (-0.3%) 

c 

Head 

No Wind 

-27 (-60) (-1.4%) 

-15 (-1.4%) 

d 

Tail 

No Wind 

+29 (+64) (+1.5%) 

+18 (+1.7%) 


These type of sensitivity studies are Important because they produce 
information necessary for the Implementation process. Computing the 
optimum profile vlth the wrong aircraft or environment models can cause 
a large percentage of the expected gain from the flight management systeni 
to not be realized. 

Many more sensitivity cases than those described above need to be 
obtained for determining the sensor and measurement processing requirements 
assoclrted with implementation of an optimum vertical flight management 
system. However, with the availability of the OPTIM and TMGEN programs, 
the user is in a position to obtain these results. 

Benefit of Free Cruise Altitude OPTIM was used to generate optimum 
profiles where the cruise altitude was both free and fixed. Figure 22 
shows a comparison of twin-jet aircraft profiles where altitude is fixed 
at 10 km (33000 ft) and free. The initial mass of the aircraft is 40.78 
tonne (90,000 lb). Range traveled is 750 n.mi. As can be seen, the climb 
and descent profiles for both cases are equivalent. The free-cruise 
altitude case begins cruise sooner (~ 32000 ft) and climbs during cruise 
to about 11.4 km (37500 ft). 
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Table 10 prcaents a cost conparison of the two profiles shown in Fig. 
22. As can be seen, the free-cruise altitude profile takes 70 sec longer 
to ccmplete, but It requires 154 kg (339 lb) less fuel than the fixed** 
cruise altitude case. (This is a reduction of 3.9J). These results were 
also based on generating optiiauiii profiles using cost of fuel of $0. 33/kg 
($0. 15/lb) and coat of time of $600. /hr. With these cost figures, the 
direct operating cost reduction of using the climb cruise is $29.20 (or 
1.37X) for this flight. 

! 

The fixed cruise results could have been Improved perhaps, by selecting 
a fixed cruise altitude of some value other than 10 km (33000 ft). How- 
ever, the pilot is not usually given that freedom. Another improvement 
could have been realized by allowing a step climb in the cruise segment 
of the fixed-cruise altitude case - e.g., from 10 km (33000 ft) to 11.28 km 
(37000 ft). An interesting and valuable addition to OPTIM would be to 
Include the step climb option. 

' Benefit of Constraining Time-of-Arrival To understand the reasons 

f why fixed time-of-arrival flight path control would be beneficial, we re- 

state three assumptions made in Chapter I concerning the future scenario 
of commercial aviation: 

1. Because of the Increasing cost and scarcity of jet fuel, air- 
craft will soon he nominally flying along minimum fuel vertical 

f flight paths. 

2. Because of Increasing demand for air travel, increasing conges- 
tion and delays of variable length will be occuring at the major 
terminal areas. 

3. Because of increasing capabilities being developed and implemented 
in communication and computer technology, the ATC system will be 
able to anticipate terminal area delay times. 11. e controller 
will be able to inform the pilot early in the flight what the 

^ expected delay will be, and he will be able to assign the pilot 

an open time slot (time-of-arrival; at the terminal feeder fix 
or outer marker. 


If these assumptions hold, the pilot will have a choice of two strategies 
to follow to take a fixed delay into account: 
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Table 10. Coac Coaiparlaon for Fixed and Free Altitude Profilea 


Range «• 

750 n.ml. 

Twin-Jet Aircraft 

W • 40.78 tonne 




° (90000 lb) 


$0. 33/kg 

- $600.00/hr. 

h ■ 10 km 

r 

($0. 15/lb) 

c 

(33000 ft) 


1 1 

Quantity 

Climb 

Cruise 

Descent 

Total 

Fuel - kg (lb) 





Fixed 

1099 (2426) 

2740 (6049) 

167 (368) 

4006 (8843) 

Free 

1091 (2409) 

2580 (5695) 

181 (399) 

3852 (8504) 

Distance - n.ml. 





Fixed 

85.6 

577.5 

86.9 

750.0 

Free 

84.6 

566.1 

99.4 

750.0 

Time - hr:m:s 





Fixed 

13. *39 

1:20:27 

15:32 

1:49:39 

Free 

13:31 

1:20:01 

17:16 

1:50:49 

Cost - $ 





Fixed 

427.73 

1530.53 

199.61 

2157.86 

Free 

424.37 

1483.67 

220.63 

2128.66 


1. Continue to fly hia nominal minimum fuel path and then enter a 

"minimum fuel flow" holding pattern to absorb the delay at the 

end of the cruise segment, or 

2. Regulate his flight path by slowing down so that he arrives at 

the terminal area within an acceptable tolerance of the assigned 

time-of-arrival . 


The algorithm developed in this study generates the optimum vertical flight 
path between a city pair which minimizes fuel and meets the delayed time- 
of-arrival constraint (Option 2 above). The fuel reduction of using this 
strategy is now compared to that of Option 1 as a function of the delay 
time . 
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Figure 22. Comparison of Free and Fixed 

Twin-jet 


Cruise Altitude Profiles for the 
Aircraft. 


Th« »«thod u.«d to . profll. that follow tht .tr«t«gy of 

Option 1 c.n b. d..crlb.d with th« .letch .hnwn in Fig. 23. It 1. ...u«.d 

t «t th« profil. follow. th« logmants b.twwn th« .wqumet of f oint. »hown 
in Fig. 23. Th« sogmnts iollow«| wroi 

1-2; Clinb along a minlmua fual sagnant. 


2-3; 


3-4: 


4-5; 


5-6; 


crulsa aagmant to tha point 
w.jara daacant would normally bagln. 

Continua at cruiaa altltuda and airspaad until tha 

mlntoum fual flow alrapaad ia obtalnad 
during tna nominal daacant. 

Dacalarata to tha mlniraxnn fual flow alrapaad whlla maln- 
aining cruiaa altltuda. rhia baglna tha holding pattern. 

the holding pattern at cruiaa altltuda and 
piriod! *irapaad to abaorb tha fixed delay time 


6-7; Continua with minimum fual daacant. 


Th. f„.l bum.d durlog „ch of th... c.n b. obt.ln.d from th. 

OPTIM progrui norm.1 printout. ThI. profU. 1. optlml.tlc in that it 
...ume. th.t th. .ircraft Xmv.. th. holding patt.rn (.t crui.. .Itltud. 
and minimum fu.l flow) with no dLcontinultp with th. optimum d..c.nt pro 



Figure 23, Sketch of Profile with Holding Pattern (Option 1), 
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Plgur«s 24 and 25 abow tha aax>unt and parcent of fuel aavad uaing 
Option 2 inatciid of Option 1. The independent variable ia the arri*'al 
tine delaye for the aediua range tri-jet aircraft with naaa of 68275 kg 
(150,000 lb) at takeoff. The range traveled ia the other variable para- 
meter in these (29 min 10 sec), about 660 kg (1450 lb) of fuel can be 
potentially saved. Approximately 225 kg (500 lb) of fuel can be saved 
for anticipated 5 minute delay, Independent of range. 

Figure 25 shows the percentage of fuel saved for the cases shown in 
Fig. 24. Up to 62 of the fuel used by Option 1 can be saved with this 
controlled time-of-arrival capability. The values shown in Fig. 25 are 
computed by dividing the reduced fuel amount by that used for controlled 
time-of-arrival (Option 2). 

Tlie results Just presented are conservative in the sense that the 
holding pattern assumed to obtain the Option 1 results is ideal. Usually, 
holding patterns are made at lower chan cruise altitudes. The results 
are optimistic in the sense that for Option 2, it 1 h assumed that the pilot 
is informed of the upcoming time delay right after takeoff. Thus, further 
study is necessary to model a more accurate representation of the holding 
pattern and to consider the cases where the pilot is informed of the delay 
somewhere during the cruise segment of flight. However, the potential 
savings are clearly Indicated. 
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Figure 25. Perceatage Fuel Saved Using Option 2 as a Function of 

Time Delav and Range. 
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IV 


SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 
Sunmary and Conclusions 

All evidence Indicates that we are rapidly consusilng our supply of 
hydrocarbon fuels. This makes It mandatory that systems be developed that 
allow conseirvatlon of the remaining sup; lies until alternate fuels can be 
developed. This Is especially true of aircraft flight whlcl' currently has 
no alternatives. This report addresses the development of on-board 
algorithms ior vertical steering of the aircraft to minimize fuel con- 
sumption and cost. 

Chapter II derives an algorithm for computing the optimum vertical 
profile using range as the Independent variable. Both fuel and time are 
penalized, and the longitudinal wind effects are taken into account. The 
Hamiltonian is constant for this mechanization, and it is equal to the 
minimum cruise cost per unit distance traveled. To obtain optimum climb 
and descent profiles involves minimizing a single function at discrete 
p-'ints along the trajectory by proper choice of thrust and airspeed. This 
algorithm proved to be a dual to the one derived by Erzberger where energy 
was used as the independent variable. A computer program (called OPTIM, 
described in Appendix A) was used to obtain the optimum vertical profiles 
for typical medium range transport aircraft flights based on these 
algorithms . 

Chapter II also presents a rethod of generating a minimum fuel flight 
path when the time-of-arrlval (or length of flight is fixed. A conver- 
gence problem sometimes occurs when using this option because aerodynamic 
and propulsion data are stored in tabular forms rather than as continuous 
functions in the programs. 

In Chapter III, the accuracy of the vertical profiles obtained from 
OPTIM are examined by using a more complete longitudinal model of the 
aircraft. This model was Incorporated into a computer program called 
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TRAGEN (described In Appendix E) which steers the aircraft to follow an 
input or computed reference trajectory. Results show that the OPTIK 
reference trajectories are both flyable and accurate in totma of fuel 
burned and time expended. 

Chapter III also pressnts the changes in characteristics of optimum 
reference trajectories due to changes in range, initial mass, and wind 
profile for a typical transport aircraft. It is concluded that these 
changes are simple modifications to the nominal reference profiles; they 
could be used to compute perturbations to a nominal profile on-board with- 
out recomputing the entire reference profile. 

Chapter 111 also examines the effects of errors in the estimated 
values of initial mass and the wind profile on the performance obtained 
during climb. A 9X increase in initial mass (6804 kg (15000 lb)) can 
cause a 23% increase in time and a 19% increase in fuel required to achieve 
the desired cruise conditions. Wind errors have a smaller effect. Tliese 
sensitivity studies are useful for specifying how accurately various para- 
meters which affect the flight performance need to be measured. 


Chapter III also illustrates further utility of the OPTIM and TRAGEN 
program capabilities. Results are given which: 

a) . Compare the fuel and time costs of a typical optimum profile 

flight where the cruise altitude is constrained to 10 km or 
is free. The range was 750 n.mi., and a twln-Jet transport 
model was used. The path with the fixed cruise altitude con- 
straint required 154 kg (3J9 lb) more fuel which is 3.9% of 
the total fuel burned. 

b) . Investigate the btnefits of having fixed t ime-of-arrival guid- 

ance capability on-board the aircraft. The fuel that could be 
saved bv using this capability to absorb delays rather than 
using holding patterns was computed for delays of up to 30 min 
and ranges of 500, 1000, and 1500 n.mi. About 225 kg (500 lb) 
of fuel could be saved for a 5 min delay of a tri-jet aircraft 
(regardless of range) using the fixed t tme-of-arr ival option. 

Fuel saved could be as high as b% of the total fuel used. 

cl. Compare handbook reference and optimum profiles. Tt was shown 
that an optimum climb for a trl-Jet aircraft would use 169 kg 
(371 lb) less fuel (7.8%) than the handbook reference climb pro- 
file. However, after making adjustments for equal range traveled, 
this savings was reduced to 52 kg (115 lb) less fuel (2.1''). 

An optimum descent requires 57 kg (124 lb) less fuel than a 
faster handbook desct>nt profile. 
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Many raora Interesting studlea can be conducted with these complementary 
computer programs. 


Re conmenda t ions 

Based on the experience obtained by review of previous work, develop^ 
ment of the OPTIM and TRAGEN program capabilities, and conversations held 
with many government and industry personnel associated with research and 
design in the a.r transportation and air traffic control fields, the 
following recommendations are made concerning immediate future work. 

The following extensions to tho OPTIM and TRAGEN programs would pro- 
vide the capability to obtain more realistic reference flight paths, 
correct known minor program errors, and add capabilities which would be 
useful for obtaining improved solutions for other phases of flight; 

a) Modeling of lift and drag 

The aerodynamic data used to model the aircraft lift and drag 
in OPTIM and TRAGEN are quite complete. However, in the 
course of the study it was found that great care had to be 
taken in obtaining lift and drag coefficients from table look- 
up routines. This was particularly true for drag; the drag co- 
efficient is currently given as a function of the lift co- 
efficient and Mach number. From a cross plot of against 
Mach number for a given a highly non-linear result occurs 
in the region of typical cruise points (e.g., M = .78 and 

= .35). Thus, tables which interpolate linearly between the 
lines of constant Mach number will create ridges or kinks in 
the data which can have an advi rse effect on the determination 
of the optimum cruise Mach number in the OPTIM program. The 
kinks produce discontinuities in attempting to converge on 
time-of-arrival . One solution would be to use curve fitting 
techniques to insure that the drag coefficient is always a con- 
tinuous function of both and Mach number. This also re- 
quires that the aerodynamic data be carefully ch >cked to 
adjust any data points causing kinks to occur. 
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b) , M odeling of propulsioD data 


Thft propulsion data used to model the turbofan engine for the 
twin-jet aircraft is strictly "uninstalled" data. For establishing 
trends and for making relative comparison, these data are suf- 
ficient. However, if one wishes to model a given aircraft more 
exactly and to predict absolute results, installation losses 
must be Included in the model. The losses include inlet pressure 
recovery losses, exhaust nozzle losses and bleed and power 
extraction. The following table gives data for one particular 
engine model to illustrate the point. 


Table 11 . Example of Unmodeled Propulsion Losses 


Mach Number * 
Altitude ■ 
Standard Day 
Reference Thrust ■ 
Reference SFC “ 


0.78 

9.1A4 km (30000 ft) 

18593 N (4180 lb) (uninstalled data) 
0.773 


Item 

Tlirust Change - = — 
N 

SFC Change - ~~ 

Nozzle Gross Thrust 
Coefficient 

-.0008 

+.0004 

Inlet Pressure Recovery 

-.0041 

+.0015 

Bleed and Power 
Extraction 

-.0462 

+.0295 

Total 

Installed Thrust « 17641 N 
(3966 lb) 

Installed Slu * 0.797 

-.0511 

+.0314 
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These qhenges are nontrivial and should be considered when modeling 
an actual system. At Idle conditions > the effects are even more 
pronounced for the bleed and power estracti'on because of the low 
level of thrust. It is recommended that these Installation effects 
be modeled and that the reLults be compared with those of the pre- 
sent study to establish quantitatively the installation effects. 

c) . Correction of the time-of -arrival convergence problem. 

Currently, because most of the aerodynamic and engine character- 
istics are obtained (for both aircraft models) from tables using 
linear interpolation, there are kinks in the data. This causes 
convergence problems when the time-of-arrival is constrained. 

Steps to correct this problem would be to 

a) . Convert aero and engine data to continuous functions 

so that interpolation can be eliminated, as dis- 
cussed above, 

b) . Ensure that there are no double minimums in the cruise 

cost at a particular altitude, 

c) . Replace or improve the optimization routine, and 

d) , Check out OPTIM with the modifications. 

This correction is vital to smooth running of the algorithm on- 
board an aircraft. 

d ) . Optimum cruise - step change in altitude . 

Currently, OPTIM is based on optimum cruise climbs. This added 
option would be an extension to using fixed cruise altitudes. 

OPTIM would determine where to change altitudes in a step fashion 
as fuel is burned off dur-*ng cruise. 

e) . Exte nd TRAGEN to simulate the entire flight path . 

Currently, TkAGEN can be used to simulate either the climb or 
the descent profile. This addition would add logic so that 
cruise (constant altitude, step climb or cruise climb) could 
be simulated along with a three-segment climb-cruise-descent 
profile . 


77 


f) . TaVLftoff parfomnanc* n>od«liag. 


Currantly, OPTIM begins l[light path optimisation calculations 
at an arbitrarily salectad spaad and altitude which is designated 
as the beginning of climb. If this program is to be used as a 
flight planning tool, some allowance must be made for the time, 
fuel burned and altitude attained during the takeoff phase. This 
phase Includes the ground roll, rotation and climb during flap 
retraction. The takeoff performance should be studied to develop 
sufficient parametric data to be able to predict the velocity 
and altitude at the beginning of climb and the fuel burned during 
the takeoff. Both the ambient temperature profile and airport 
altitude must be taken Into account. 

g) . iiescent performance modeling. 

Currently, OPTIM computes the descent portion of the trajectory 
with the engines usually at a flight idle power setting. In doing 
so, two potential constraints are Ignored - the rate of Increase 
in cabin pressure and the maximum allowable fuselage angle. Both 
factors are Important for a commercial transport aircraft. The 
descent trajectory for both aircraft should be studied to de- 
termine if these constraints are significant in relation to the 
overall flight path optimization methodology used in OPTIM. If 
these constraints do prove to be significant, it is necessary 
to include these effects in both OPTIM and TRAGEN. In addition, 
the fuel, time, range, and altitude covered during the final des- 
cent phase with flaps and landing gear deployed need to be de- 
termired. These effects would also be included in the OPTIM 
and TRAGEN programs. 

h) . Cruise control to account for wine and temperature 

Currently, OPTIM uses a single wind profile which specifies 
magnitude and direction as a function of altitude. This profile 
is used to compute the optimum climb, cruise, and descent por- 
tions of flight. There is an eventual need to develop variable 
profile models of the wind and temperature as functions of both 
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altitudtt and horizontal poaitlon« Than, by using thaaa models, 
the optimization techniques need to be modified to account for 
the modeled wind and temperature variations. This is especially 
important for the cruise phase of flight which may cover a 
large range over which the wind changes may be considerable. 
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APPENDIX A 


TRAJECTORY OPTIMIZATION USING THE ENERGY STATE METHOD 

Th« purpose of this Appendix is to summarise briefly the theore- 
tical background used for optimisation of trajectories. These principles 
form the basis of the numerical process used for computing the optimum 
vertical profile of a Jet aircraft. More details are given in Refs. 17, 

21, and 23. Reference 25 derives the principles upon which trajectory 
optimisation is based. In Refs. 17 and 21, Ersberger and Lee apply these 
principles using the energy state approximation to obtain a practical, 
efficient means of generating the optimum vertical profile. 

OPTIM is an extension of the original computer code developed by 
Ersberger and Lee that encompasses these principles and is based on their 
methods. The latter part of this appendix describes how OPTIM is organized. 

In the following sections, the theory of trajectory optimization is 
first presented. Then, the application of this theory to minimizing the 
direct operating cost (DOC) of an aircraft traveling over a fixed range 
is outlined. This is followed by a discussion of the details of going 
from theoretical expressions to a practical computer code. The theore- 
tical points are presented without proof, for conciseness. The reader 
wanting more detail should review the references. 

Theoretical Principles 

In Ref. 25,, a description is given of the requirements for solving 
an optimization problem in\olving a continuous dynamic system with nc 
terminal constraints but with fixed terminal time. Ihis description is 
repeated here because it presents the basic principles which extend to 
the aircraft profile optimization problem. 


81 


W!:ce/>1NG fioT 


A systma (th« aircraft) ia govamad by tha nonlinaar differential 
equations 


X • f(x,u,t) ; 


x(t^) given; 

^ i ^f * 


(A.l) 


where x is tha n-dimansional state vector and u is the m-dimensional con- 
trol vector. Tlia cost function which is to be minimized is of tho form 


<|>(x(tf),tf) + y* ^L(x,u,t) dt. 


(A. 2) 


Here, (]) Is the terminal cost function, and L is the cost per unit time along 
the trajectory. The problem is to find the sequence of controls u(t) that 
minimize J. 


First, the system equations are adjoined to J with the multiplier 
vector A(t): 


J • <t>(x(t^) ,t^) + 


y ^L(x,u, 


t)+X (t) {f (x,u,t)-x) } dt. (A. 3) 


Then the Hamiltonian function is defined as 


H(x,u,t) ■ L(x,u,t) + X (t)f(x,u,t). 


(A.4) 


Equation (A. 3) is integrated by parts to yield 

J » <{>(x(t^) ,t^) - X^(tj) x(t^) + x"^(t^) x(t^) 

+ y * {H(x,u,l) + X^(t) x(t)} dt. 


(A. 5) 


Next, the change in J due to variations in u(t) and x(t) is con- 
sidered for fixed t and tr.‘ 

o f 


1 


6J 


/ + )«* + IS«“ j <“• 


U.6) 


Th« •l«m«ntf of X(t) arc chosan to caus« tha coafflclanta of £x In Eq. 
(A. 6) to vaniah undar tha intagral and at t^i 


»T M It 3f 

" ■ 3x * " 3x ■ 3x ’ 


- If • 


(A. 7) 


Equations (a. 7) are called the co-state aquations. Than, Eq. (A. 6) becoaas 


i 5J ■ X^6x(t ) + 
o 


/' 


{ 3H 
3u 


6u dt. 


(A.8) 


For J to be minimum, 6j must be zero for arbitrary u(f); this implies that 
for no bounds on u, 


in 

3u 


» 1 t < t^ 


on the optimum path. If the control variables are constrained as 


n(u,t) V 0, 


(A. 9) 


(A. 10) 


then for u(t) to be minimizing, we must have <5J 0 for all admissible u(t). 

This implies, from Eq- (A.8) that 


P iu CH > 0, 


(A. 11) 


for all t and all admissible 5u(t). In other words, H must be .minimized 
over the set of all possible u; this is known as the minimum principle [25]. 


In summary, to solve for u(t) that minimizes J, the differential 
equations (A.l) and (A. 7) must be solved simultaneously, where u(t) is 
determined from Eqs. (A. 9) or (A. 11). The boundary conditions on the state 
X at t^ and > at t^ are specified, resulting in a two-point boundary-value 
problem. 
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If L and f ar« not axpllclt function! of t, than 


H 


3h . 3H • 

^ ^ X + 


3H . 3H • 

^ u + 


3h 




u + 




(A. 12) 


Each clamant of Eq. (A. 12) li zero on the optimum trajectory, from \dilch 
we can concliide that H la constant on the optimum trajectory. This latter 
point is used in the analysis presented in Chapter IX. 


Application to Aircraft irofile Optimization Uning the Energy 

State Approximation 


Here we are concerned with applying the above theory to the pro- 
blem of choosing the thrust and airspeed values to control the aircraft 
vertical profile in going from one point to another. The cost function 
J is the direct operating cost (DOC) 'hich is the sura of fuel and time 
costs. This is, in Integral form, 


J 



w + Cj.)dt 



(A. 13) 


wtieft- ( ^ is the cost of (uel (S/kg or $/lb), w is fuel flow (kg/lb/nr), C^. 
is the cost of time ($/hr), i> the direct operating cost, and t^ is the 
time to fly the specified distance traveled d^. It is also assumed that the 
typical vertical profile is as shown in Fig. A.l - that is, it contains 
climb, cruise, and descent portions which have the constraint that 



d , 1 

dn 


(A. 14) 


where 

d * x(t .) ■ the distance traveled from the start print to 
up cl ‘ 

where t le cruise segment begins (at time t * t ,). 

c 1 

d, * d, - x(t -) “ the distance traveled from the end of 
dn f cf 

cruise (at time t ■ t^^) to where the descent segment ends. 
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Tim*, th« co.t function (Eq. (A. 13)) 


can be rewritten a* 


■/ c^dt + w, - ^ y« 

t-, 

/■ 




(A.15) 


+ (d,-x(t^^)- id^-x(t^^]), + yt, 

where 1 « the cost per unit distance while In cruise 
Simplified point 


C^dt. 


cf 


craft are 


•"'“‘‘“o- of longitudinal motion of 


the alr- 


*■ (T-0)/m - g ^ ^ 

h - VgSln y , 

X ■* V cos Y + V 

a y » 


(A. 16) 


where the flight path angle (y) dynamics and mass loss d 

are neglected. Here, ^ hum 


'^a * nlrspeed, 
h » altitude. 


- longitudinal oomponont (tangant to earth'a aurfaca) of 
m - aircraft mass. ^urtace) of 


wind speed. 


T * thrust, 
D * drag. 
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Al.sc, tt Is assumed that lift L * mg cos Tne eltect of ituias loss is 
accounted for by cont inucus.ly updating mass without adding another state 
variable. 

Specific energy K is defined a& 

E * b + vJ'/Zf’ , U.17) 

£t 

whicl. is the sum of potential and hineVlr eriergy per liolt mass. Its 
time derivative is. fcand to be 

£ » V (r"D)/r.ig . (A. 18) 

•A 

Tlit; enorg ' state appioximation is bused on the assumption that, potential 
and kiiU'i -c energy can be interohan'j.od instantaneously. In this approxima- 
tion, Che energy state variable replaces altlttide and airspeed state 

* A 

variables U-b 1 . Thus, Eq. (.A. 17) can be used in place of and h in Eq, 

(A. 1 6) . 

It is assumed that the aircraft specific energy Increases ir.onotonically 
during i.'iimb end decreaaef.? nionotonically during descent. Tliis assumotion 
ill u.sea in the development to ciiange ti’e i ndependent variable in Eq. (A. 1.5) 
trom time to ertergy. This use.s the transformation 

gt = ^ IA.19) 

L 


It is mathematicallv a.onven i one to evaluate the lact integral in Eq. 

(A. 15) backwards in time so that the energy state Is monotonically increas- 
ing during its evaluation. This means that the running distance (range) 
variable during the descent can be measured backwards from the end point. 
Thus, we can think of range measured in two ways as shown in Fig. A. 2. 
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CRUfSF 


I 




Figuic A. 2. Measurcnu^nt oi Range from the Oripln or fo the Destination. 


In this sketch. 


’'up^'cP' 
»dn<'> ■ 


range measured on the way up in forwtird time t, 
value of when initial cruise Is reached, 
range measured on the way up in backward time x, 
value of when final cruise is reached 


Also, we define the variable, x to be range traveled during climb and 
descent. The di stance traveled during cruise is then constrained to be 
(d^ ~ x)- t'<e can tiuni sec that an incremeatal change dx in the range 

variable x is eouivalent to inf'n-.ment-al chances in both x and x. . 

“ up dn 

That is 


dx 


d(x 


up 


* ’‘dn> 


(A. 2 


From this discussion, the second of Eqs. (A. 15) can be written as 


J -/''iCjdt + (df - * +y^ct IcJdT . 


(A. 2 


We use Eq. (A. 19) and the transformation 



to rewrite Eq. (A. 15) as 




J>0 




dll' cf 


/‘•'(fl 


CA.23) 


Here, ’ ^'ct ’ valoec of “ju’rp’y ffaltmted 

at tltne t equal to 0 Oiici a.ul sMiif’ r evaluated at v ^ r.nd • . rcspw-'C” 


rlvely. 


Note £looi Etj . (A. 23) that t.hf rjrij'.e variable x only appujjs ar fho 
sum of clltub and dosreut distances •*- Thus, the state equation 

for this system oi equations can be vrrltten as 


(i (x + X , ) 
dx ^ up dll 

iiE dE” 


/fV 4 V I \ / (V 4 V ) \ 

. { I . ( l„dn wd>./ \ 

\ E /e>o I lyf A 


(A. 2^0 


here, and are the .lonijitud hiai cemponents oi the v/imi sf.eed 

for climb aad desc *nt. rhea, analogous to Eq, f,',.'!', tlit i?anil,tcinlrai ii 



1 \ 


'(V -fV )\ /( 

V , ; t \ 1 ~ 

y 

) " 

\ i \ ) ’ ) ' 

np yup ^ 

wdij' j j 

\ F 

I ft 

E>0 

'E<0 ' 

C ' 

E>0 

'e<o' 


(A.25) 


This can be divided as 


“ ■ 


■1 

C . + X (V +V ) 
d up wup 

A 


Pl» 

1 

•T 

E>0 

1 

•rzi 

1 


(A. 26) 


J E<0 


Now, analogous to Eq. (A. 7), the costaCe equation for ,\ can be written 


as 
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ilfO 




cA 


ax 


an 


d(x +x, ) 
' up dn 


« 0 


(A. 27) 


and from Eqs, (A. 7) and (A. 23), this costate has the final value 


. 3(ld-~x 

3(> f up dn 


CA.28) 


where 'j' is the cruise cost per unit distance. 


Note, this problem could be placed in a slightly more conventional form 
by dividing it into two problems - one for climb and one-half of the 
cruise distance and the other for decent .in<( rhc other half of the cruise 
distance. Then Eqs. (A. 27) und (A. 28) would be replaced by 


’3E 


IH 

3>; 


(A.25) 


up 




5([d,/2 - x^l - - 


dx 


’P 


up 


for climb. i''or doscend, j 

i 

(A. 30) \ 


This allows splitting the Hamiltonian defined in Eq. (A. 26) and allows for 
ii)(E^^) ^ fact, in the actual implementation ^ E^^ because 

optimum cruise energy changes as fuel is burned off. The principal results 
are unchanged , however . 

Thus, from Eq. (A. 11), (A. 29) and (A. 30) the trajectory optimization 
problem becomes 


:r\ dll 

. — H- ja. ^ 

dE d:< . 

dn 


0 


Dad^/2 - 


X(Ecf) 


cf^ 


3x 


dn 


- 'P (E,f) 
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CA.Sl) 



Thus, the optimization problem reduces to rolviag pointwise mlnliBum 
values of the algebraic functions dellnod by Bq. (A. 31) during the climb 
and descent portions of the traji; 'story. 

Equations (A. 29) and (A. 30) are #oa transve reality condition for the 

free final state problem (d c, < d-) with terminal cost. Thus, the 

up dn r 

constant value of A for climb ar- 1 descent is found to be the negative of 
the cost per unit distance for crulae. 

The cruise cost ij; (« - A) If found by assuming that the aircraft is 
in static equilibrium during cvui-'^e (T ■ D), and that 


■MEc) 


mlu __d 

V (V +‘~v"T 

C C w 


(A. 32) 


In other words, for any cruise altitude, there is an optimum thrust and 

airspeed such that the cost per unit distance minimized. The 

optimum cruise cost as a function of cruise energy is typically of the 

shape shown in Fig. A. 3. Thus, there is also an optimum cruise energy 

Ecopt where cruise cost is minimized. If the range is long 

enough so that there is sufficient range to reach optimum cruise energy 

E it should be done, and the cruise conditions should be set so that 

copt 

• ♦<®copP- 
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Figure A. 3. Optimum Cruise Cost as a Function of Cruise Energy 


For the case where there is no cruise segment (d^. = 
the cost function contains only integral terras. Then, the transversaJ Ity 
condition yields V «■ - if‘(t^). That is, X would be the negative of 
ijj(t^), where is the optimum tost for cruising at the highest point 

reached on the climb trajectory. 


The optimum cruise energy is only specifically reached when 

there is range enough to climb to and descend from the optimum altitude/ 

airspeed values, where C (E^) is minimum. For ranges less than this 

value, the maximum value of F that is reached is a free variable less than 

c 

the optimum value. its choice is made to optimize the cost function of 
Kq. (A. 2 3). 


From Eqs. (A. 23) and (A. 25), one can write 


3J 

3E 


H + 



^dn> 


>*(E^)] 


3e' 


0 , 


at E * E . This is 
c 


(A. 33) 


91 


t 


(A. 34) 


H + d 
c 


ay; 


- 0 


where le the cruise distance , and is the total value of H (H^p + 

at the cruise point. Thus, Bq. (A. 34) can be used, along with other 

characteristics of ij/ and H, to determine the relationship between E^, and 

d . The Hamiltonian evaluated at E * E is the cost penalty to achieve a 
c c 

unit increase in cruise energy. For 0, Eq. (A. 34) can be written as 


“c 

Figure A. 4 shows the family of trajectories which have this characteristic. 

These occur at values of E below E . where 3<{'/3E < 0 (see Fig. A. 3). 

c copt 

That is, non-zero cruise segments occur at short ranges with cruise energies 
less than the optimum energy value for long range. 


For the case where H •• 0, d is zero for 3iJ;/?E < 0. The distance 

c c 

can be non-zero only at optimum cruise energy where 3ij^/5E “ 0. This 
family of trajectories is shown in Fig. A. 5. 


Thus, we have a situation where positive values of dictate one 
type of trajectory and zero valuet dictate another. In Ref. 17, it is 
shown that if the aircraft engine specific fuel consumption S„„ is in- 
dependent of the thrust T (so that w ■ S_„T) , then the structure of the 

r L> 

trajectories will be like Fig. A. 5 with no cruise segment occurlng except 
at (This J.mplies that the Hamiltonian is zero at the maximum energy 

point). For this case, the optimum thrust setting for climb is T^^^, and 
the optimum setting for descent is 

If the engine specific fuel consumption is dependent on thrust, and 
the thrust values are not constrained during climb or descent, it is 
shown in Ref. 17, that the Hamiltonian H^ is again zero at the cruise 
energy, and again the trajectory structure is like those of Fig. A. 5. 
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Figure A. 4. Optimum Profile Energy ve Range for > 0 at Cruise. 



Figure A. 5. Optimum Profile Energy vs Range for H » 0 at Cruise. 



If tha Is depsiulsnt cm thrust, sod coastrslnsd to ths maximum 
vslus for climb sod to ths minimum idlt vslus for dsscsnt, th«n ths 
Hsmlltonisn is positivs st cruiss. This csusss positivs cruiss ssgmsnts 
according to Eq. (A. 35) at cruiss snsrgiss bslow ths optianmi. For this 
cass, ths optimum trajsctoriss will havs shapss similar to Fig. A. 4. 

Thsss trajsctoriss ars slightly Isss sfficisnt than thoss of Fig. A. 5. 
bscauss ons Isss control is available for optimiaaticm. 

Some Mechanization Details of ths Computer Program 

The remaining ssctions of this Appendix describe iiow the previous 
theoretical material has been utilized to construct an offline computer 
program for generating optimur vertical profiles for models of medium range 
tri-jet and twin-jet transport aircraft. This material is presented In 
an alternate way in Ref. 21, and the program is referred to here as 
OPTIM. 

By examining the specific fuel consumption data of the turbojet 
engine, it is determined that is dependent on thrust. Thus, for the 
transport models, two types of short range profiles must be considered - 
those represented by Fig. A. 4 (.Type 1 profile) when thrust is constrained 
and airspeed is the single control - and those represented by Fig. A. 5 
(Type 2 profile) when both thrust and airspeed are used as controls. 

The solution to optimum climb and descent profiles is found by 
minimizing the Hamiltonian expressed in Eqs. (A. 31). The independent 
variable (energy) is stepped along in fixed increments (e.g., 150 ra 
(500 ft)), and the Hamiltonian is minimized at each energy setting. 
Minimization occurs by finding the best values of airspeed 
(and optionally thrust climb function and the 

descent function are individually minimized. 

To solve Eqs. (A. 31) requires knowing two more quantities: 

X or >[i(E ) the cruise cost per unit distance. This comes from 
^ evaluating Eq. (A. 32) at the desired cruise altitude. 
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- th« crulsa «n«rgy. Thi« 1« a function of th« crula* altitude 
and tha aaaoclatad crulaa alrapaad obtained in Eq. (A. 32). 

Mota that for tha Type 2 profile at abort rangea, there ia no cruise 
segment. In this case* the maximum energy achieved at maximum altitude 
is referred to as the cruise energy E^. At that altitude, there still 
is defined a minimum cruise cost according to Eq. (A. 32). 

For the Type 1 trajectory of short range, there exists a non-sero 
cruise segment which is determined by use of Eq. (A. 35). To solve Eq. 

(A. 35) requires that the Hamiltonian defined by Pqs. (A. 31) be solved 
at the point of transition from climb-to-cruise. It also requires 
knowing the slope 3i{;/3E of the cruise cost for a change in craise energy at 
that point. 

Cruise Optimization 

The first step that must be taken to compute optimum trajectories 
is to derive the optimum cruise cost and its derivative 3ip/DE. This 
is done by computing what is referred to as the "cruise table". The 
parameters that affect this table are the assumed cruise mass, the 

t 

wind profile, and the lift L, drag D, thrust T, and fuel flow w 
characteristics of the aircraft. The optimization process searches over 
the acceptable ranges of altitude and airspeed for a given mass. The 
results are collected in tabular form for a series of different assumed 
cruise masses. 


Again, the minimum cost of flight during cruise per unit distance 

for a fixed cruise mass W is found by 

c 


X 




min 

V 

a 



w + 



+ V 


w 


(A. 36) 
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This Msunas that tha aircraft ia in static aquilibrlini during crulsa, i.a.. 


T cos a - D, 

L + T sin a « ng, 


(A. 37) 


whara tha angla-of-attack a is found by solving thasa aquations sinul- 
tanaously. Tha altituda is stappad in 305m (1000 ft) Incraaants from saa 
laval to calling altituda (whara maximum thrust Just balancas drag) . At 
altltudas balow calling altituda, tha alrspaad ~ dapandant drag curva 
crossas tha maximum thrust curva at two points (V^ and Vj) as lllustratad 
in Fig. A. 6. Thus, for aach altituda laval, tha ^ ,luas of Vj^ and V 2 
ara datarmlnad, and than ^(W^,E^) is minlmlsad with raspact to alrspaad 
batwaan thasa two limits. Kastrlctlons ara that ba graator than 

0.1 Mach and that ba less than 0,89 Mach for tha trl'-jet (0,84 for the 
twin Jet) for structural reasons. 


After the cruise cost is minimised at each discrete altitude level, 
these numbers are stored in a table with altitude as the independent 
variable. Typical results .,>'*e plotted in Fig. A. 7, Presented here are 
also the optimum cruise Mach number M and the optimum thrust setting 

opt 

EPRopj.. After results are obtained in steps of 305 ra (1000 ft), the 
minimum cost point is found as a function of altitude. In the OPTIM pro- 
gram, the cruise table optimization results are obtained by using a 
Fibonacci search with ten Fibonacci numbers. (See Ref. 22). 


The cruise table results are obtained for cruise mass varying as 
dictated by the program input. Usually, the cruise mass is incremented 
in steps of 2268 kg (5000 lb). Up to ten values of cruise mass can be 
used. For each cruise mass, the optimal cruise altitude, cost, speed, 
power setting, fuel flow and specific energy are computed. An example 
of optimum cruise cost and a function of cruise mass is shown in Fig. A. 8. 

C limb Optimization 

After the cruise tables are generated, the program proceeds with 
obtaining the optimum climb trajectory. This requires guessing what the 
cruise mass will be, based on the takeoff mass. The guess is used to 
obtain a trial value for (or X) in the Hamiltonian from the cruise 
tables. The procedure to obtain this guess is based on an empirical for- 
mula which iterates until convergence is made. 
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F igure A.7 . 


Optimum Cruise Cost as a vunction of Altitude for 
Mass of 54432 kg (120,000 lb). 
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Cruise 
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Figure A. 8. Optimum Cruiaa Cost and Cruise Altitude as Functions of 
Cruise Mass for a Tri-jet Aircraft Flying Into a par- 
ticular Head Wind. 


The climb optimization process starts by assuming (E^) “ 1. 5 4' , 
where is first obtained by setting the initial cruise mass W . 

COpt JO 

equal to the takeoff mass (an input). The appropriate cruise tables 
are used to int^’rpolate to find the corresponding value of E^ associated 
with 1.5 Then, an empirical formula of the form 

''up ■ “-I' 

is used to obtain an approximation to the fuel burned to reach E^. Here, 

E^ is the takeoff aircraft energy, is a reference mass (61690 kg 

(136000 lb) for the tri-jet) and is the previous value of cruise mass. 
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Then, the cruise mess is updated at - F^p. This process is 

repeated until the difference in consecutive estimates of P falls below 

up 

45 kg (100 lb). 

When the cruise mass estisiate is obtained, the corresponding values 

of E and i>(E ) are obtained from the cruise tables. Then, the program 
o c 

is ready to generate points on the optimum climb trajectory. This is done 
by incrementing specific energy and minimising the Hamiltonian function 



(E) 


V ■»• C^. - (V^ + V^ ) 


(A.39) 


at each point. (This is the first of Eqs. (A. 31)). That is, the program 

2 

starts with initial energy * h^ /2g. It steps the energy a fixed 

amount AE (,say 150 m (500 ft). At this point, it searches over true air- 
speed V (and possibly thrust setting tt) so that Eq. (A.39) is minimized. 

8L 

For the turbojet engines, thrust is governed by EPR settings which vary 
between 1.1 (idle thrust) and some maximum value leas than 2.4. The true 
airspeed has an upper limit governed by 


a) . 0.89 Mach structural limits (tri-jet); 0.84 (twin jet) 

b) . 250 kt (IAS) below 3048 m (10000 ft) for ATC restrictions, 

c) . V 2g (E-h) which insures that the aircraft climbs, and 

d) . V^, the upper value shown in Fig. A. 6 where max thrust 

equals drag. 


Tlie lower limit is governed by 

a) . Vj^, the lower value shown in Fig. A. 6 where max thrust 

equals drag, 

b) . 0.1 Mach 

c) . 1.5 m/sec (5 ft/sec) less than the previous value of V 

to limit large jumps in flight path angle. 

The Fibonacci search technique is again used to determine V and tt 
which minimize Eq. (A.39) for the fixed value of energy E. The value 
chosen for airspeed is accurate to within .0056 Mach, and EPR is accurate 
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to within .009. Associated with these values of and it are values of 
energy rate E (Eq. (A. 16)) and altitude h: 

h - E - V.^/2g . (A.AO) 

From these, approximate values of time, range, flight path angle, and 
fuel burned are obtained from 

At - AE/E, (A.41) 

slny ■ (Ah/At) /V , 

cL 

X ■ EAx ; Ax “ (V cos Y + V (h)) At 

^ A W 

F - ZAf ; AF - w At. 

The above process Is repeated by stepping along energy In Incre- 
ments of Ae until Is reached. The last value of Eq. (A. 39) Is 
stored for possible use In evaluating the cruise distance. 

The above climb optimization procedure Is repeated with 'J' ■ 1.5'P^, 
l.Oijj^, and perhaps other values until the total range of flight converges 
to the appropriate value. This Is discussed in further detail later. 

Descent Optimization 

The descent optimization is very similar to the climb optimization 
with regard to the equations which are evaluated. The optimization pro- 
cess requiies estimated values of and at the beginning of descent, 
and an estimate of mass at the end of descent. The method used 
to obtain these estimates is discussed in the next section. 

If there is a cruise portion of flight, fuel will be burned during 
cruise. Thus, the value of and at the beginning of descent 

will usually be different than at the beginning of cruise. If there is no 
cruise portion, then these values will be Identical. 

The descent profile is obtained by starting at the final energy 
state and then going backwards in time. The energy rate is constrained 
to be negative with respect to forward time. 
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similar descent profile constraints exist on true airspeed as for 
those of the cUmb profile* The thrust level Is on or near the Idle value 
during descent. 

Cruise Fuel Bum 


To estimate the final mass during cruise (U^^) and landing (W^), 
the following steps are taken: 


1) . Determine 4*^1 the Initial cruise cost based on the Initial 

cruise mass obtained from the climb optimization. 

2) Use the Initial cruise mass and t'^ compute the fuel flow 

3) . Estimate the cruise range by the empirical equations, 

P - - 1.5 . (A.42) 

c copt 

d ■ b, + b„P^ + b»?^ b,P + be . 

A). Compute the cruise fuel as 

d^/(V + V„(h^)). (A. A3) 

c c c c w c 


5). Estimate the average cruise mass as 

W - W - 0.5F 
c c c 


(A.AA) 


6) . Use the cruise Cable Co obtain the corresponding cruise cost 

e __ 

altitude h^, fuel flow w(ip^), true airspeed V^, and wind 
speed V^(h). 

7) . Recompute Eq. (A. A3), and then find the final cruise mass. 


W 


cf 



F . 
c 


(A.A5) 


8) . Use the value W , in the cruise tables to obtain ip(W ,). 

cf ct 

As with the climb, set ip « 1.5 >jj(W^j). 

9) . Use this value of to obtain h^^ and E^^ from the cruise 

tables. These are the end conditions for the descent trajectory 
obtained backwards in time. 
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10). Estimate the landing mass from tha ampirlcal formula 


P - 1.5, (A. 46) 

Wj - W,., - (c^P^ + CjP + Cj). 


The values of 4*, E^^ and obtained by the above procedure are used for 
obtaining the optimum descent t’'...ti>- ory. The descent portion of the 
Hamiltonian Is of the form 


C, w + C. 


“dn<« 


- ♦«cf> 


(V. + 


V ) 

W 


(A.47) 


this function Is also minimised at each of the given values of energy. 


After the first descent profile Is completed, a new estimate of 
cruise distance is obtained by using Eq. (A. 35), or 

^^c “ . (A. 48) 

Then, step (4) above Is repeated to obtain an Improved cruise fu&l burn. 
Then, the improved landing mass estimate Is 


W. 


W. - (F 
i 'up 


+ F + F . ) 
c dn 


(A. 49) 


The landing trajectory is reoptimized with this new value of 
landing mass. Then, improved values of total range traveled, time 
required, and fuel burned during climb, cruise, and descent are made. 


For short range flight, the above steps ass'.uned that a Type 1 trajec- 
tory is generated because thrust is constrained to maximum value during 
climb and idle value during descent. If thrust is free, then a Type 2 
trajectory will result, with no cruise portion. For this case, the steps 
required to estimate cruise distance d^ and final cruise cost, mass, and 
energy cap be eliminated. 
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Cruise Coat Estimatlnn 

Tb. flr.t dirt, .nd d..c.nt proflU. .r. E«..r.cd with - 1.5* 

(E ). Ih* next net le generated with - l-O* <'copt'‘ ° 

th“e» «lue. of cmlee co.t heve «. eeeocUted tenge on the curve ehown 

in Fig. h.9. M th* “tel tenge de.lted 1. greeter then (the velue 

^ 4 1 Oib (E )). then it le aeeumed that the optimum cruise 

obtained using leOy '“copt^^* . . 

” t. j •riaanm a third set of climb and descent 

altitude end energy are reached. Then a third aet cii. 

/_ V e-hi B rase the cruise distance 

profiles is generated using 'l'(Ecopt^* ^ ^ ud j 

is computed so that the desired overall range is exactly achieved. 


If the desired range is between end Rj^^^ 

a ■ \ *>teA A 


in Fig. A. 9, then an 

Iterative procee. 1. need to ohteln end the a.eoclated de.lred 

ra.ge. Iteration, are .topped when the totel r«.ge traveled 1. within 
.one a.«U dl.tance t of the de.lred range. (In OPIIM. t 1. ..t at 5 n. 


mi. ) 


103 


Cost Parameter Ratio - (^ / v(E 





Figure A. 9. Relationship between the Cruise Cost Parameter and 
the Associated Range of Flight 
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AFPEilDrX B 


ENGINE MUPEL DEVELOPMENT 


The modeling of the turbofan engine for the twln-Jet aircraft in 
program OPTIM required changes to subroutines FCLIMB, ENGEPR, and ENGIDL 
[22]. Idle thrust and idle fuel flow are computed in ENGIDL and FCLIMB, 
and thrust and fuel flow at all other power settings are computed with a 
call to subroutine ENGEPR. The actual computations are done in subroutine 
ENC.EP3. All engine data is stored in tables in the BLOCK DATA sub- 
program . 


Idle Performance 


Fipure B.l shows idle performance (net thrust and fuel flow) for 
the engine taken from the engine installation handbook. Note that surge 
bleeds can be either opun or closed depending upon altitude and flight 
Mach number. To model this phenomena, two sets of tables were developed - 
one for bleeds closed and one for bleeds open (see Figs. B.2 and B.3). 
Then, an additional table was developed to determine the altitude for 
surge bleed closure as a function of flight Mach number (see Fig. B.4). 
Figures B.5 and B.6 show idle performance for cold and hot days. It is 
assumed that the fuel flow varies in inverse proportion to the square 
root of the T.j .2 ratio. This is expressed as 


'^fuel non STD 
, 

"fuel STD 


1 ^12 non STD 

J T^2 


(B. 


1 ^ 2 ’ degrees absolute, is the stagnation temperature at the compressor 
inlet. There is essentially no change in idle thrust due to temperature 
variation. 


105 





o 



l-l 

x: 

I 

• » 

I 

o 

lu 


D 

U. 



Altituc 

Figure B,2 Engine Idle Performance 
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ICAO Standard Atmoaphara 
Surga Blaada Cloaad 


Min Fuel Flow ■ 245 kg/ hr 
(540 Ib/hr) 





Fuel Flow - w, Ib/hr 


M - 1.0 


ICAO Standard Atmoaphare 
Surge Blacda Open 




Altitude - km 


Figure B.3 Engine Idle Performance - Surge Bleeds Open 
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Figure li.-4 
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Ho Sorvie* Hl««d Mo {*ow«r Extraction 

Standard Atnoapharlc Conditlona - 40*F XOOZ Ea« Eacovary 

ICAO Standard Ataosi^ara 



Surxc Blaadi Open 



Altitude - km 


Figure B.6, Estimated Engine Performance at Flight Idle on a Hot 

(+22*K) Day 
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KftxiaiuM Etiiia* Pr«»tur« latlo 


HaxiNun •ngltM prMsur* ratio ^EPt) is dtlarminsd at aacb flight con- 
dition to insurt that tha angina is not sat bayond its oparational liaits. 

As shown in Fig. B.7« tha naxiatum EPU for climb is a function of tha con- 
prassor inlat taaparatura (T^2^ (which is cosiputad from Mach number 

and altituda). Likawisa, maximum EPR at cruiaa is a function of 7^2 only. 
Two diffarant curvas ara used: ona for altituda £ 9.144 km (30000 ft) and 

ona for altituda > 10.668 km (35000 ft). Bar/aan 9.144 and 10.668 km 
(30000 and 35000 ft)» and tha maximum EPB is assumad to vary linaarly. 

Engine Thrust 

Figure B.8 gives corrected thrust (thrust/amblant pressure ratio) 
as a function of flight Mach number and tha engine pressure ratio for the 
twln-jat aircraft. Figure B.9 is a repeat of Fig. B.8 with additional 
dat^ which repreacntn the same engine performance in the tri-jet aircraft. 

It can be seen that with one exception the tri-Jat data lies below the 
twin-jet data.* It was found that the tri-jet engine data represented 
installed thrust with the values being an average of the three engines on 
the aircraft. Because the middle engine on the tri-jet aircraft suffers 
greater Installed losses because of the "S" shaped inlet duct, it is not 
surprising that the thrust shown in less. The program selects the appro- 
priate engine data for each aircraft. 

Engine Fuel Flow 

Engine corrected fuel flow is shown in Fig. B.IO as a function of the 
engine pressure ratio and altitude. Note that the corrected fuel flow is 
also a function of flight Mach number for EPR less than 1.6. The correction 


* The one exception xS at M ■ 0.7 and EPR ■ 1.85. This data point was 
presumed to be in error, and a smooth curve is assumed. 
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lure Ratio, tPR 



Corrected Thrust 










Corrected Thrust 









Corrected Fuel Flow • w/*.^ * 10 /hr 


i 


I 

c 

« 

u. 



Figure B.IO, Corrected Fuel Flow of Engine 
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factor for fu«l flow, K^, ■ aipplo liuaar equation developed from the 

straight line shown on the inset in Pig, B.IO. That is, 

- .0022 (T^2^ +0.97 . 

with T^2 degrees centigrade. 
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APPBiroiX c 


MODELLIKG OF THI TWTS-JET AIRCRAFT ARROPTHAMIC8 


Program OPTIM ccmputas th« dynamic forces in flight for tha complete 
trajectory - climb, cruise and descent. Thus, aerodynamic models are re- 
quired to compute lift and drag forces over the complete range of altitude 
and Mach number in the flight envelope of the aircraft. The following 
sections outline the twin-jet aerodynamic models now used in OPTIM. 


Drag Force 


The dimensionless aerodynamic drag force coefficient is separated 
into its important contributing elements. At a given angle of attack, a, 


“ SbASIC ^^DGEAR 


(C.l) 


where is the basic drag coefficient for the airplane in free air, 
with the landing gear retracted, no spoilers deflected and not In ground 
effect. The curves include trim data for level flight. 


The low speed drag coefficient is shown in Fig. C.l as a function of 
angle-of-attack and flap position. The flaps-up data on this plot are for 
stall only, where C, is a function of angle-of-attack. The flaps-up data 
are shown only for angles-of-attack larger than 8 degrees. 

The high speed drag coefficient, showing Mach effects for the flaps-up 
configuration, is plotted in Fig. C.2 as a function of Mach number and lift 
coefficient C^* (C^^* is defined later) with gear retracted, no spoilers de- 
flected and out of ground effect. The M £ .6 curve is only to be used for 
lift coefficients smaller than .757, which agrees with angle-of-attack of 
8 degrees or less for a trimmed airplane. 
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FPCCCOJHG BLAfM tiCT HllVtt'.r 


Basic Drag - 




Basic Drag with No Flaps 



Id Eq, (C.l) tb« cocfficirat lncr«gMsit In bMlc drag 

coafficlant du« to th« main and noaa landing gaar axtonaion givan by tha 
following: 

^^DGEAR " ^^^SxJKAR^MK) 

Tha low-apaad drag Incraraant, ^ Pl-ottad in Fig, C«3 x > a 

lunction of angla-of-attack and flap aatting, Tha flapa*nip Mach nunbar 

factor, ^^DGEAR^M la plotted in Fig, C.A. 

^^DGEAR^>W) 

Specific data points are identified :^n each of Figs. C.l - C.4 des- 
cribed above. These points are tabulated in subroutine CDRAC3 in Program 
OPTIM. 


Lifr. Force 

The dimensionless aerodynamic lift force coefficient of the airplane 
is separated into its important contributing elements as 


C 


"BASIC 


r&i 1 

L\d.Ti /m " \ dft /m-0 J • 


(C.3) 


Also, 




(C.4) 


GEAR 


Here, C, * Is used in the calculation of the basic drag coefficient. C 

c ^'BASIC 

is the basic lift coefficient for the rigid airplane at low speed in free 
air, and with landing gear retracted. The coefficient is plotted as a 
function of angle-of-attack ot, and flap setting in Fig. C.5. Also, the 
term Rc ) - (C ) I is the deviation of the basic airplane flaps-up 

L S M 0 M-Oj 
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lift csMtficlmt lit a • 0* with Kwh m»b«r. Th« t«r« 
is plotttd oo Plf, C«6. 


Ths tsrsi» [(SI -(SUl a Is ths dsvlstlon from ths bssic lift 
cosfficisnt dus to ths vsristion of ths lift curvs slops st a ■ 0* with 

Msch Dimbsr snd sltituds. Ths tsns [(ttIm -(a^)»w)] is plottsd in 


Fig. C.7, snd ths >H) vslus is ths slops of ths flsps up 



St a - 0*. 


As with ths drsg dsts in ths prsvious ssction, spscific dsts points 
srs idsntifisd in Fig, C.5 - C.7, Thsss points srs tsbulstsd in subroutins 
CLXFT3 in Progrsm OPTIM, 
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''nSASIC'M 



IncremeutJ*’ Change in Lift Curve Slope 



Altltud* - km 



Figure C.7 Incremental Cliangcs xn Lift Curve Slope as a Function of 

Mach Number and Altitude 
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i 

1 
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APPIWDIX P 


CLIMB FUEL B8TIMATB 

Fuel bum«d during climb is ftstlmatcd In subroutine FULEST [22] so 
that the aircraft mass at the beginning of cruise can be determined before 
actually computing the climb portion of the flight. A new correlation of 
the climb fuel flow was established using data from OPTIM runs. 

Figure D.l is a plot of the correlation which shows that the fuel 
fraction of initial mass used in climb is a function of the change in 
energy state ~ ^co^ initial gross mss. 

This correlation is a? follows: 

<=cl - 

(D.2) 

(D.3) 

numerical values, Che OPTIM program uses 

« fuel burned in climb, lb, 

> Initial climb mass, lb, 

» true airspeed, ft/sec , 

» altitude, ft , 

- 32.2 ft/sec^ . 
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Here, for these 

”fCL 

W 

to 


a 


h 

g 


SLOPE - 2.9 X 10“^^(Wj.^) + .3625 x 10“^ , 


where 


2g 


4* h • 







APPESDIS E 


AIRCRAFT EQUATIONS OP MOTION AND AUTOPILOT MODET.S 


The objective of the TKAGEN prograw Is to eiinulate a transport air- 
craft being steered to fly along either an input or computed reference 
trajectory. This trajectory may either be a climb or descent profile. 

The simulation must be accurate enough such that the performance of the 
aircraft (in terms of fuel burned anJ time required to reach the destina- 
tion point) is adequately determined, as measured from the output. 

Adequate accuracy is obtained with a five-state variable longitudinal 
aircraft model. 

The purpose of this appendix is to present the analytical expressions 
upon which the simulation was developed; this is done in two parts. The 
first section below defines the overall system and presents the differential 
equations of motion and fuel bum. The second section describes different 
methods for generating typical guidance commands and autopilot equations. 


Equations of Motion and Fuel Burn 

To examine the vertical pt-oflle of the aircraft (i.e., altitude and 

airspeed vs range), the longitudinal equations of motion are of primary 

importance. The short period equations of motion and the throttle dynamics 

are Ignored. Thus, the control variables in this longitudinal plane are 

the angle-of-attack a and the magnitude of the thrust vector T. These 

quantities are shown with respect to aircraft true airspeed V , lift L, 

a 

drag D, mass W, and flight path angle 7 in Fig. E.l. 

The kinematic equations of motion of the aircraft in the longiLivdlnal 
plane are 

X “ ^ '^w * (E'l) 

h >» V sin y , 

a ’ 
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Figure E.l Vector Diagram of Longitudinal Fbrces 


where 


X - distance, or range, measured on the ground, 
h - altitude, 

^ wind speed. 


The inertial speed along the true airspeed vector V is 

A 

- V + V COSY 
law 

From Fig. E.l, the time rate of change of this vector for constant Y is 

• • • 1 » 

“ V + V cos Y “ ~ (T cos a - D - W sin y)« 

1 a w ' ra ' 


The time rate of change of the wind speed is 


V 

w 


Substitutiug Eq. 
V 

a 


3V 

Dh ^ ’ 

V sin Y 
3h a 

(E.4) into Eq. (E.3) and solving for 

3V 

— (T cos a >- D - W sin y) ~ V 

m dh a 


V leaves 
a 

sin Y cos Y 


(E.2) 


(E.3) 


(E.4) 


(E.5) 
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Also, from Fig. E.l, on* can writ* 

K - - (L CO* Y - W + T Bin ( Y + a ) - D *in Y ). (E.6) 

81 

Equation* (E.5) and (E.6) repraaant tha fcinatic aquation* of motion of 
tha aircraft. 

Tha ramaining tarm that mu*t ba accountad for i* tha tima-varying 

maa* of tha aircraft. Spacifying tha thruat alao apacifia* tha fual 
• • 

flow w. Thua, tha maaa changaa at tha rata 

W - -w . (E.7) 

Equations (E.l), (E.5), (E.6), and (E.7) are the five basic equations 

used to represent the longitudinal dynamics of the aircraft. Lift and 

drag (L and D) are computed as functions of V , h, and a . Fuel flow 

8 

w is a function of V , h, and thrust T. 

a 

Further refinement could be added to these equations to include the 
effects of the following: 

1) . throttle dynamics (including transient fuel flow); 

2) . relationship between throttle position, EPR setting, and thrust; 

3) . short period dynamics relating time rate of change of angle- 

of-attack, pitch rate, and pitch angle to elevator deflection; 

4) . required turning (lateral) motion for flying over fixed waypoints; 

and 

5) . lateral wind and gust effects. 

However, these effects are considered to be of second order, and not re- 
quired for the intent of this simulation. For a more exact autopilot 
simulation, they would be required. 


133 


Th* flight path angla la dafinad aa 


Y - ain"^ (E.8) 

By dlffarantating thla axpraaalon and ualng Bqa. (B.5) and (E.6)i ona obtalna 

1 2 

Y ■ (T ain a - W coa y + !■ + m ain y), (E.9) 

Equation (E.9) can ba uaad In placa of Eq. (E.6). 


Staerlng Procedurea 

The reference trajectories which are given to be followed consist of 
a sequence of points containing values of time, range, altitude, airspeed, 
flight path angle, specific energy, mass, and other variables. Any of 
these quantities which is measurable and monotonically changing can serve 
as the independent variable. To minimize airborne computer memory re- 
quirements, It is important to make the stored data representing the re- 
ference trajectory as compact as possible. 

In this study, a set of steering equations is used to take points 
from the reference trajectory, convert these points to reference trajectory 
commands, and then use these commands to set values of the control vari- 
ables. This steering process represents a rudimentary form of an autopilot. 

The steering process consistc of commanding the thrust T and angle- 
of-attack a values so that the aircraft follows the reference as closely 
as possible. The system that Includes this process is depicted by the 
block diagram in Fig. E.2. Note that flying along a reference trajectory 
consists of steering to connect a series of reference points. When a re- 
ference point is reached, new steering commands must be Issued so that the 
aircraft will then be guided to the next reference point. 


134 


1 



Figure E>2. Elenents of the Longitudinal Aircraft Model 


To fly along the reference path* an independent variable ie first 
chosen, For this study, two different independent variables were chosen - 
altitude for climb and range for descent. Then, the remaining variables - 
primarily airspeed, flight path angle, and altitude (for range as the in- 
dependent variable) - are stored as tabular functions of the chosen in- 
dependent variable. 

Also, it is possible to fly along a nominal path using two approaches: 

1) . An open-loop approach where the thrust vector is directed 

in such a way over the next period that by the end of that 
period the next reference point is reached. 

2) . A closed-loop approach where the aircraft is continually 

steered to a continuously commanded trajectory which connects 
the reference points. 

Both of these approaches were examined for simulation of flying the climb 
profile. However, good results were only obtained with the closed-loop 
approach. 
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A probltm with op«n-loop •tMtlng is that It aaaunta that conataat 
or litiaarly varying controls will causa tha and points of a rafaranca 
prof 11a to ha connactad. This assumption doas nut account for parturba- 
tlons dus to wind, ate. along tha way. Although tha opan-loop nathods 
produca paths which hava roughly eorract valuas of alrspa«d aiul altltuda 
at glvan rang* valuas, thara vara larga axcurslons from tha iTSfaranca 
flight path angla for the climb profiles. 

Another problem with tha open-loop approaches was that both a and T 
ware varied to achieve fixed values of and h for given range points. 

For optimum climb, thrust Is usually set at the maximum value. Thus, 
usuall" only a remains as a valid control variable. 

Another consideration for Implementing the climb profile is that 
there is no reason why a particular cruise condition (altitude, airspeed) 
has to be achieved when a certain range x Is reached. Thus, a more logical 
Independent variable is altitude, with range allowed to be a free variable. 


For these reasons, a closed-loop steering approach was devised where 

reference values of flight path angle (with respect to the air mass) and 

airspeed are obtained as functions of altitude. (This assumes that 

altitude is monotonically Increasing during climb.) A perturbation control 

law was set up so that variations In a from a reference value a were 

o 

proportional to variations in y and V from their respective command values. 

d 


Because y and V tend to change linearly with time, they can be 
considered as ramp functions. Thus, the closed-loop controller should be 
considered to be at least a Type 1 system. From Eqs. (E.5) and (E.9), with 
no wind, the system perturbation equations are 


mSO 

a 

m V (Sy 

a 


T sin a 6a - 


5a 


T cos a 6a + 


12 

9a 



6a 



6Va - W cos Y<5 y» 

6V + W sin y5y* 
a ' 


(E.IO) 
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Th« rasultlng traiuifar fuoctloiui batwaan Y» 



Gi » •*• ^ 

(a/oj)^ + 2; (a/w) + 1 


iot 


(a/u)^ + 2C(a/u) + 1 


whara the ti»a cooatanta and othar parainatara are functiona of tha para 
matara in Eq. (E.IO). 

Tha control problam can now ba intarpratad aa ahown in Fig. B.3. 

To obtain the Type 1 ayatam, tha control law haa to ba of tha form 

K ^4 

to-(Kj+-^) - V.) + (Itj + 



V 

nom 


Y 

nom 


Figure E. 3 Control Loops for Perturbation Control of 
Airspeed and Flight Path Angle. 


( Ell) 


(E.12) 
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wh.r. and art tha cow^ndad valuaa of V and Y. Thia la tha i < i 

Proportlonal-plua-intaaral controll.r r 4 * ^ ^ claaalcal 

da* i rad chotan to produca tha 

d«.irad ra.ponaa for r-oval of proflla arrora. 

To ganarata tha contlnuoua conaumda V, and Y durlna cH«h m, 
tiona »ada at aach rafaranca point ara ’ co«puu- 


P+1 

^n+1 " 


V - V 

a n+1 ''ar 

iT. , - h 

n+1 n 


'n ^ . 


'^a “ V + (h - h ) . 

c *n P 3h • 


Wlien the flight path angle la very small (durlna th. 1 

fHah^ u taurxng the initial period of 

flight and when the aircraft level* off at 3048m (10000 ^ 

fore resuming climb), Eqb (k 14 > do ► i ° speed be- 

rao;, tqs. (t.i4) do not work well. For these cases 
more appropriate to set cases, it is 


a » 

^ n +1 


and use the control law 


Ja . (Kj + -1) (y_ . 


(E.16) 


J 


Tb« above mchod worked quite well in ceueing tho simleted profile to 
closely follow the reference path. Only one set ol' gain values was 
sufficient for the entire trajectory. 

For decanding flight, the thrust again is usually constraint (idle) 
for optinun performance. Also, for this case, the main concern is to reach 
a fixed altitude when range»to~go to the destination point is a certain 
value. Thus, above 3048 m (10000 ft), the airspeed can be allowed to be a 
free variable. For this case, only inertial flight path ie required to 
be controlled. 

To generate a continuous commend Yj,,. the computation made at each 
reference point is 


Yi^ - - V ^ <Vl - V ] ■ 

Then, the control law is similar to Eq. (E.16), l.e. 

8 ^ ^'^IC “ 

where inertial values of flight path angle ar** used rather than those with 
respect to the air mass. Equations (E.17) and (E.18) form the basis for 
closed- loop control of descending flight. Again, one set of gains is 
sufficient for the entire descent profile. 


139 


WEfKRIlKIBS 


1. Covay, R.Ii» Macc«tti» GtJ.t iu)d RosMlsr* W.U., **lxMiiuitlon of Co»- 
Mrcial Aviation Enargy Conaarvatlon Stratagian/' Raport No. ATR-79 
(776D-1, Aarospaca Corpi, El 8agundo» CA. Octobar 1978. 

2. Banka. C.H.. "Ovarviaw of Tachnology AdvancaaMnta for Bnargy Efflclant 
Tranaporta," AIAA Atnoapharlc Flight Machanlca Confaranca. Papar Ho. 
79-1651, Bouldar, CO. Auguat 1979. 

3. Stangal, E.P., and Marcua, F.J., "Bnargy Managanant Tachnlquaa for 
Fual Conaarvatlon In Military Tranaport Alrr.raft," AFPDL-TR-75-156, 
Analytic Sciancaa Corp., Raadlng, MA, Fabruiry 1976. 

4. Aggarwal, R., at al., "An Analyala of Fual Conaarvlng Oparatlonal 
Procaduraa and Daaign Modificationa for Bonbar/Tranaport Aircraft," 
AFFDL-TR-78-96, Vol. II, Dynamlca Raaaarch Corp, Wllolngton, MA, 

July 1978. 

5. Yarrington, A.L., and Burghart, J.T., "Tamlnal Conflgurad Vahlcla 
Fual Conaarvatlon Sanaltlvlty Study," Boalng Commarclal Airplana Co., 
Doc. No. D6-42939, Saattla, WA, July 1976. 

6. Collins, B.P., "Eatimatlon of Fual Conaunptlon of Commarclal Jet Air- 
craft from Path Profile Data," SAE International Air Tranaportatlon 
Meeting, Cincinnati, OH, May 1980. 

7. Smith, B.A., "Automated Fuel Managment Set for Teat," Aviation Week 
& Space Technology , April 21, 1980. pp. 191-193. 

8. Klaas, P.J., "f.ontroller Numbers Keyed to Advances," Aviation Week & 
Space Technology , April 28, 1980, pp. 61-63. 

9. Sorensen, J.A. Morello, S.A., and Erzberger, H., "Application of 
Trajectory Optimization Principles to Minimize Aircraft Operating 
Costs," Procedings of the 18th IEEE Conference on Decision and Con- 
trol, Ft. Lauderdale, FL, December 1979, pp 415-421. 

10. Sorensen, J.A., "Concepts for Generating Optimum Vertical Flight Pro- 
files," NASA CR 159181, September 1979. 

11. Schultz, R.L., and Zagalsky, N.R., "Aircraft Performance Optimization," 
Journal of Aircraft , Vol. 9, No. 2, January 1972, pp. 108-114. 

12. Erwin, R.L., Jr., et al., "Strategic Control Algorithm Development," 
Boeing Commercial Airplane Co,, Contract No. DOT-TSC-538, December 1973, 

13. Knox, C.E., and Cannon, D.G., "Development and Flight Test Results of 
a Flight Management Algorithm for Fuel Conservative Descents in a Time 
Based Metering Traffic Environment," 18th IEEE Conference on Decision 
and Control, Ft. Lauderdale, FL, December 1979. 


141 




^WCEDINQ PACE BLANK NOT TOMEO 


14. TobU«, h., Q'Biri«n, P,J,» «q 4 faliMT, E.A., ''SjjiuXAtioa Study of tho 
Effoct of Pual-Consorvatlv* Approaches on ATC Procaduras and Tsnsinal 
Araa Capacity,” 8AK Air Tranaportatlon Maatlng, Boaton, H4» May 1978. 

13. Eutotmkl, E.8., "Enargy Approach to tba CaiMiral Aircraft Parfonsaiica 
Problan,” J. of Aaronautleal gciancaa. Vol. 21 Mo. 3, Iterch 1954 
pp. 187-19TI 

16. Bryaon, A.E., Daaai, M.N., atad Hoffiaan, W.C., "Enargy-Btata Approxl** 
mation of Suparaonlc Aircraft,” .1 . of Mrcraft. Vol. 6, Mo. 6, Mov- 
Dac 1969, pp. 481-488, 

17. Erabargar, H., and Laa, H.O., "Charactariatica of Conatrainad Optlnaigi 
Trajactorlas with Spaclfiad Ranga,” RASA TM 78519, Sapttndjar 1978. 

18. Erabargar, H., McLaan, J.D., and Barman, J.R., "Flxad-Ranga Optlmin 
Trajactorlas for Short Haul Aircraft,” RASA TR D-8115, Dac«abar 1975. 

19. Barman, J.B., and Erabargar, H. , "Plxad-Ranga Optimum Trajactorlas 
for Short-Haul Aircraft,” J. of Aircraft , Vol. 13, Ro. 10, Octobar 
1976, pp. 748-754. 

20. Shoaaa, H., and Bryaon, A.B., "Alrplana Minimum Fual Flight Paths for 
Fixed Ranga,” SUDAAR Mo. 499, Jtanford University, Stanford CA, March 
1976. 

21. Laa, H.Q., and Erabargar, H., "Algorithm for Fixed Ranga Optimal 
Trajactorlas,” NASA Technical Paper 1565, 1979. 

22. Anon., "OPTIM - Computer Program to Generate a Vertical Profile which 
Minimizes the Aircraft Direct Operating Costs,” Revision No. 2, MASl- 
15497, Analytical Mechanica Assoc., Inc., Mountain View, CA, October 
1980. 

23. Anon., "TRAGEN - Computer Program to Simulate an Aircraft Steered to 
Follow a Specified Vertical Profile,” Revision No. 2, NASl-15497, 
Analytical Mechanics Assoc., Inc., Mountain View, CA, October 1980. 

24. Coykendall, R.E., Curry, J.K., Domkc, A.E., and Madsen, S.E., "Study 
of Cost/Benefit Tradeoffs for Reducing the Energy Consumption of the 
Commercial Air Transportation System,” NASA CR-137891, June 1976. 

25. Bryson, A.E., and Ho, Y.C., "A pplied Optimal Control ,” Blaisdell, 
Waltham, MA, 1969. 


